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Abstract 

This  thesis  demonstrates  the  design  of  a  lateral  flight  controller 
for  the  0135  aircraft  using  the  frequency  response  approach  of  Profes¬ 
sor  Izaac  Horowitz,  University  of  Colorado,  Boulder,  Colorado.  A  single 
controller  is  designed  for  the  entire  flight  envelope  of  the  aircraft 
from  high  speed,  high  altitude  cruise  to  a  low  altitude  approach  condi¬ 
tion. 

A  three  degree-of- freedom  lateral  model  is  reduced  to  a  two  degree- 
of-freedom  model.  The  equations  of  motion  assume  a  linearized  aircraft 
model  operating  around  equilibrium  flight  conditions  allowing  only  small 
perturbations  from  equllibruim. 

Robustness  of  the  design  is  guaranteed,  since  the  variation  in 
flight  parameters  are  considered  throughout  the  entire  design  process. 
Robustness  is  demonstrated  by  simulating  the  compensated  aircraft  at  the 
extremes  of  the  flight  envelope  considered. 

This  thesis  concludes  that  the  approach  of  Professor  Horowitz  is 
useful  for  the  design  of  flight  controllers  for  the  C-135  aircraft. 
Some  software  development  to  aid  in  the  design  process  is  recommended. 


MULTIPLE  INPUT  -  MULTIPLE  OUTPUT  FLIGHT  CONTROL 


DESIGN  WITH  HIGHLY  UNCERTAIN  PARAMETERS 
APPLICATION  TO  THE  C-135  AIRCRAFT 


I.  Introduction 


Dr.  Isaac  Horowitz  of  the  University  of  Colorado,  Boulder,  Colorado 
has  developed  a  synthesis  technique  for  multiple  input  -  multiple  output 
feedback  control  systems.  He  has  worked  with  the  U.S.  Air  Force  on  a 
number  of  projects  such  as  the  F-16  CCV.  His  frequency  response 
approach  to  controller  design  has  produced  successful  designs  for  a 
number  of  difficult  flight  control  problems  (Ref  9:5-8). 

The  Air  Force  is  interested  in  the  use  of  Dr.  Horowitz'  approach 
for  a  variety  of  aircraft  control  problems.  This  thesis  deals  with  the 
design  of  a  lateral  flight  controller  for  the  C-135  aircraft  using  his 
technique . 

This  first  chapter  presents  some  background,  a  statement  of  the 
problem,  the  scope,  assumptions,  approach,  and  the  sequence  of  presenta¬ 
tion  for  the  rest  of  the  thesis. 

Background 

An  aircraft  in  flight  is  an  uncertain,  non-linear,  highly  coupled 
multiple  input  -  multiple  output  (MIMO)  control  system.  The  engine(s) 
and  the  flight  control  surfaces,  elevator,  ailerons,  rudder,  flaps, 
spoilers,  and  often  others  provide  the  inputs  to  the  system.  The  out¬ 
puts  are  the  velocity  vector,  the  angles  and  the  angular  rates  of  the 
aircraft  body  with  respect  to  some  predetermined  axis  system.  In 


general,  each  input  exists  to  primarily  affect  only  one  or  two  of  the 
outputs.  But  in  practice,  an  input  usually  affects  a  number  of  outputs 
in  addition  to  those  it  was  designed  to  control.  For  instance,  the 
ailerons  control  roll  rate  and  thereby  control  bank  angle.  But  in 
addition  to  a  roll  rate,  an  aileron  input  produces  a  yaw  rate  and  a 
pitch  angle  change  due  to  asymmetric  drag  and  a  change  in  the  lift 
vector  respectively.  Thus  the  inputs  and  outputs  are  Interrelated  or 
coupled.  These  relationships  of  Inputs  to  outputs  are  not  linear  and 
are  seldom  .known  with  a  high  degree  of  certainty  (Ref  9:1). 

Classically,  engineers  assumed  their  mathematical  models  of  these 
input/output  relationships  to  be  linear  and  invariant  in  small  regions 
around  equlllibrlum  points  (Ref  4:154).  The  MIMO  problem  was  then 
uncoupled  to  achieve  a  set  of  single  input  -  single  output  (SISO) 
problems  (Ref  3:162).  The  designs  resulting  from  this  approach  were 
often  successful  but  usually  involved  a  great  deal  of  trial  and  error 
(Ref  9:1). 

More  recently,  control  designs  consisting  of  feedback  gains  have 
been  derived  from  the  state  space  representation  of  the  aircraft  model 
using  the  digital  computer  to  perform  the  complex  matrix  operations 
required.  These  techniques  are  often  grouped  under  the  heading  of 
Modern  or  Optimal  Control.  (Ref  3:523-62).  The  robustness  in  these 
designs  tends  to  break  down  when  required  state  information  is  unavail¬ 
able  or  has  to  be  estimated.  Another  approach,  based  upon  the  state- 
variable  approach  is  the  method  of  Professor  Brian  Porter.  This  method 
has  demonstrated  ability  to  achieve  robustness,  but  can't  predict  the 
range  of  uncertainty  over  which  the  design  will  be  valid  in  advance. 

The  essence  of  Dr.  Horowitz'  method  is  the  consideration  of 
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the  performance  specifications  and  the  uncertainty  in  the  plant  quanti¬ 
tatively  from  the  very  beginning  of  the  design  process.  The  uncertain 
MIMO  plant  model  can  be  derived  directly  from  the  linear  differential 
equations  or  from  the  state  space  representation  of  the  MIMO  system.  By 
the  use  of  elementary  matrix  operations  the  MIMO  system  is  separated 
into  SISO  systems.  The  designer  derives  a  controller  for  each  of  the 
resulting  SISO  loops  based  on  the  design  specifications  and  any  loop 
design(s)  already  completed.  If  each  individual  SISO  design  meets  its 
requirements,  then  the  combination  of  them  is  guaranteed  to  meet  the 
overall  design  specifications  (Ref  11:977). 

Problem 

The  purpose  of  this  thesis  is  to  apply  the  method  of  Dr.  Horowitz 
to  a  lateral  control  design  for  a  large  transport  aircraft,  specifically 
the  C-13S.  The  primary  objective  is  that  an  engineer  with  a  reasonable 
backgound  in  control  theory  and  in  the  frequency  response  method  of 
analysis  should  be  able  to  follow  the  design  steps  and  be  able  to  apply 
this  technique  to  similar  problems. 

Subordinate  objectives  are  as  follows: 

1.  Achievement  of  competence  and  insight  in  feedback  theory, 
l.e.,  the  means  and  cost  of  attaining  desired  performance  despite  uncer¬ 
tainty. 

2.  Analysis  of  design  in  achieving  robustness,  or  the  range 
of  flight  conditions  over  which  the  design  is  valid. 

3.  Qualitative  comparison  of  the  design  resulting  from  this 
thesis  to  the  design  resulting  from  a  concurrent  0135  thesis  using  the 
method  of  Professor  Brian  Porter. 
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This  thesis  is  limited  to  a  lateral  flight  control  design  for  the 
0135  aircraft  using  the  Horowitz  method.  Comparisons  are  made  between 
these  results  and  those  of  a  Porter  method  design  for  the  same  aircraft, 
but  the  comparisons  are  qualitative  only. 

The  analysis  of  design  robustness  is  limited  by  the  availability  of 
valid  data  representing  the  flight  characteristics  of  the  aircraft  under 
varying  conditions. 

Assumptions 

The  following  assumptions  are  made: 

1.  Flat,  non-rotating  Earth  Inertial  reference  frame. 

2.  Quasi-steady  airflow  (Ref  2:22). 

3.  Rigid  body  aircraft. 

4.  Constant  aircraft  mass. 

In  addition,  it  is  assumed  that  computer  simulations  in  the  com¬ 
puter  aided  design  program,  TOTAL,  provide  realistic  representations  of 
aircraft  motion  in  flight. 

Approach 

The  first  step  in  the  approach  to  this  problem  is  to  reduce  the 
aircraft  data  into  a  matrix  representation  of  the  plant  model  at  each 
flight  condition.  The  result  is  a  plant  matrix,  J?,  such  that  £  -  lb 
where  jr  is  the  output  vector  and  £  is  the  input  vector.  Inverting  the  £ 
matrix  produces  the  effective  SISO  plant  transfer  functions  used  during 


each  loop  design. 

The  second  step  is  to  derive  a  frequency  domain  representation  of 


the  desired  specifications.  This  is  admittedly  the  weakest  part  of  this 
thesis.  Military  specifications  for  transport  aircraft  are  not  very 
demanding.  Therefore,  in  this  thesis,  much  stricter  demands  are  imposed 
on  the  aircraft  response.  These  demands  are  somewhat  arbitrary  and  may 
be  unrealistic. 

The  third  step  is  to  apply  the  Horowitz  method  to  the  design 
Itself.  The  primary  objective  of  this  thesis  is  the  presentation  of 
this  application  of  this  method.  The  method  also  promises  to  provide 
insight  into  the  costs  and  trade-offs  Involved  in  a  design  problem.  As 
these  become  apparent,  they  are  identified. 

When  the  design  is  complete,  the  fourth  step  is  to  simulate  the 
compensated  system.  This  is  done  using  the  program,  TOTAL.  Aircraft 
response  and  control  surface  movements  to  a  step  command  for  a  30-degree 
bank  turn  are  presented. 

The  fifth  and  final  step  is  to  compare,  qualitatively,  this  design 
with  the  design  resulting  from  a  concurrent  C-135  thesis  using  the 
Porter  method.  The  ease  of  application,  success  in  meeting  specifica¬ 
tions,  and  the  resulting  robustness  of  the  two  design  methods  are  the 
primary  comparisons  made. 

Presentation 

This  thesis  consists  of  seven  chapters.  Chapter  II  presents  the 
single  input  -  single  output  (SISO)  design  theory  of  Professor  Horowitz. 
The  design  of  multiple  input  -  multiple  output  (MIMO)  systems  is  covered 
in  Chapter  III.  Chapters  II  and  III  are  co-written  by  Lt.  Brian  Pawlow- 
ski,  GE-83D,  author  of  another  thesis  using  Professor  Horowitz'  tech¬ 
nique  and  this  author.  Chapter  IV  contains  the  C-135  aircraft  model 


used  in  this  thesis.  Two  designs  are  presented  in  Chapter  V.  The  first 
design  is  presented  in  detail,  the  second  is  simply  shown  for  com¬ 
parison.  The  results  of  simulations  for  both  designs  are  found  in 
Chapter  VI.  Finally,  Chapter  VII  consists  of  conclusions,  recommenda¬ 
tions,  and  a  comparison  of  the  results  of  this  thesis  to  the  results  of 
a  concurrent  C-13S  thesis  using  Porter's  method. 
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II 


Introduction 

Chapters  II  and  III  present  an  overview  of  the  Quantitative  Feed¬ 
back  Synthesis  Technique  used  in  the  design  of  multiple  input  -  multiple 
output  flight  control  systems  for  this  thesis.  Examples  are  presented 
to  aid  in  the  understanding  of  the  material.  The  technique  is  valid  for 
the  general  n-by-n  case.  However,  for  simplicity,  the  examples  below 
are  either  single  loop  or  two-by-two  systems.  For  a  discussion  of  the 
three-by-three  case  and  extrapolation  to  the  more  general  case,  see 
Reference  8  in  the  bibliography. 

The  flight  control  problem  involves  a  multiple  input  -  multiple 
output  (MIMO)  plant  requiring  regulation  and  control  due  to  parameter 
uncertainty  and  disturbances.  The  mathematical  equations  describing  the 
motion  of  an  aircraft  are  highly  non-linear.  For  design  purposes,  these 
equations  are  linearized  about  a  point  in  the  flight  envelope,  or  flight 
condition.  Uncertainty  arises  as  the  linearized  coefficients  vary  with 
airspeed  and  altitude  (see  Chapter  IV). 

The  Quantitative  Feedback  Synthesis  Technique  developed  by  Dr. 
Isaac  Horowitz  uses  feedback  to  achieve  closed-loop  system  response 
within  performance  tolerances  despite  plant  uncertainty.  The  range  of 
plant  uncertainty  and  the  output  performance  specifications  are  quanti¬ 
tative  parameters  in  the  design  process  (Ref  8:81).  The  fundamentals  of 
the  design  method  are  presented  in  the  discussion  of  the  single  input  - 
single  output  design  problem  of  Chapter  II.  The  multiple  input  - 
multiple  output  design  procedure  is  described  in  Chapter  III,  using  the 
fundamentals  developed  in  Chapter  II. 
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Problem  Definition 


The  general  single  Input  -  single  output  (SISO)  problem  involves  a 
plant  transfer  function,  P,  with  uncertain  parameters  (gain,  poles,  and 
zeros)  known  only  to  be  members  of  finite  sets.  The  design  specifica¬ 
tions  dictate  the  desired  response  of  the  plant  to  inputs  and/or  distur¬ 
bances.  The  problem  is  to  obtain  a  controller  forcing  the  plant  output 
to  satisfy  performance  tolerances  over  the  range  of  plant  uncertainty. 


The  basic  SISO  control  loop  structure  is  shown  in  Fig  II-l. 


Fig  II-l  Two  Degree-of -Freedom  Control  Loop 


In  this  figure,  r(t)  is  the  command  input  to  the  system  and  d(t)  is 
a  disturbance  input  to  be  attenuated.  P  is  the  plant  transfer  function, 
the  characteristics  of  which  are  not  precisely  known.  The  compensator, 
G,  and  the  prefilter,  F,  are  to  be  designed  to  force  the  system  output, 
y(t),  to  be  an  element  of  a  set  of  acceptable  responses,  despite  the 
uncertainty  in  P  and  the  disturbance  input,  d(t).  The  plant  input 
signal,  x(t),  is  identified  since  it  is  generally  of  interest  because  of 
physical  or  practical  constraints.  The  signals,  r(t)  and  y(t)  are 
assumed  measurable  quantities  and  the  latter  is  available  for  feedback. 
Access  to  both  signals  allows  the  use  of  the  two  degree-of-freedom 
structure  of  Fig  II-l  and  provides  the  designer  with  two  independent 


compensator  elements,  F  and  G  (Ref  9:13).  It  Is  also  assumed  that  r(t), 
y(t),  and  (for  now)  P,  such  chat  y(t)  *  Px(t),  are  all  Laplace  trans¬ 
formable  functions  (Ref  9:8). 


There  are  four  transfer  functions  of  interest  in  Fig  II-l,  where 
the  loop  transmission,  L,  is  defined  as;  L  •  G  P.  The  system  output  due 
to  the  command  and  disturbance  Inputs,  respectively,  are; 


\  -  Y(s)/R(s)  -  F  G  P 

dl  +GP) 


F  L 

ttttt 


(11-1) 


T_  -  Y(s)/D(s)  -  1 

U  (1  +  G  P) 


1 

(1  +  L)  (II-2) 


and  the  plant  input  due  to  the  command  and  disturbance  inputs,  respec¬ 
tively,  are: 


**  -  x(s)/r(s)  -  n 

(1  +  C  P) 

h  ’  W/DM  *  TT-S™ 


F  G 

7TVTT  (n-3) 

-G 

(1  +  l!)  (II-4) 


The  design  specifications  may  impose  constraints  on  any  or  all  of  these 
transfer  functions,  but  for  the  purpose  of  this  example,  only  the  first 
two  are  considered. 


Design  Specifications 

The  design  specifications,  or  closed-loop  system  response  toler¬ 
ances,  describe  the  upper  and  lower  limits  for  acceptable  output 
response  to  a  desired  input  or  disturbance.  Any  output  response  between 
the  two  bounds  is  assumed  acceptable.  The  response  specifications  must 


be  determined  prior  to  applying  the  design  method.  Typically,  response 
specifications  are  given  in  the  time  domain,  such  as  the  figures  of 
merit  M  ,  t  ,  t  ,  and  K  based  upon  a  step  forcing  function,  (Ref  3:346) 

p  8  p  m 

or  as  a  bounded  region  as  shown  in  Fig  II-2.  Response  to  a  step  input 
is  a  good  initial  test  of  system  response.  Bounds  (Tt )  and  (T  )  of  the 
figure  are  the  acceptable  lover  and  upper  limits  of  a  system's  tracking 
performance  to  a  step  input.  Desired  system  response  to  a  step  disturb¬ 
ance  generally  requires  maintaining  the  output  below  a  given  value,  thus 
only  an  upper  bound  is  necessary  as  shown  by  curve  (TQ)  in  Fig  II-2. 
Additional  similar  bounds  are  needed  if  other  inputs  are  to  be  con¬ 
sidered. 


Fig  11-2  Time  Domain  Step  Response  Specifications 


The  design  technique  is  a  frequency  domain  approach,  therefore  the 
time  domain  specifications  must  be  translated  to  bounds  in  the  frequency 


domain.  Desired  control  ratios,  ■  [Y/R]._  and  Tur.  “  [Y/D]._,  are 

MK  MK  HU  HU 

modeled  to  satisfy  the  performance  specifications  using  the  pole-zero 
placement  method  as  described  in  Section  12-2  of  Ref  3.  For  response  to 
a  step  input,  a  third  order  model  with  one  zero  is  suggested. 


A  (s  +  zY) 


M 


(s2  +  2£a>s  +a>_2)(s  +  pO 
n  n  J 


(11-5) 


The  pole-zero  pattern  of  Equation  II-5  is  shown  in  Fig  II-3.  The  loca¬ 
tions  of  the  roots  are  adjusted  until  the  step  response  of  the  modeled 
control  ratio  matches  the  bound. 


X 

jfc> 

I-  plane 

X 

Fig  II-3  Third-Order  Control  Ratio  Pole-Zero  Pattern 


The  frequency  domain  characteristics  are  considered  during  the 
response  modeling.  It  is  desireable  to  keep  the  magnitude  difference 
(as  a  function  of  frequency)  between  the  upper  and  lower  bound  models  of 


iTMR(j<u)  |  as  large  as  possible  at  all  frequencies.  Choosing  a  lower 
bound  model  with  a  greater  pole  to  zero  ratio  than  the  upper  one  ensures 
that  the  magnitude  difference  approaches  Infinity  in  the  limit  as  o 
approaches  Infinity. 

Errors  made  during  this  modeling  process  manifest  themselves  in  one 
of  two  ways.  First,  if  the  worst  acceptable  response  model  is  not 
really  acceptable,  the  system  may  not  meet  the  specifications  over  the 
assumed  range  of  uncertainty  in  P.  And  second,  if  the  entire  range  of 
allowable  outputs  is  not  considered,  the  bandwidth  of  the  compensation 
will  be  higher  than  necessary,  increasing  the  cost  of  the  compensator 
(Ref  9:5). 

Once  control  ratios  are  obtained  for  each  time  response  bound,  a 
magnitude  plot  of  the  frequency  response  (Bode  plot)  for  each  T^jto)  is 
made  on  the  same  graph  as  shown  in  Fig  I 1-4.  These  plots  are  the 
frequency  domain  representation  of  the  design  specifications  on  and 
Tq.  These  derived  frequency  domain  specifications  are  used  to  obtain 
the  bounds  on  the  loop  transmission,  L(jcu),  as  described  later. 


Nichols  Chart 


The  primary  tool  used  In  the  design  of  the  compensator  elements,  G 
and  F  is  the  Nichols  chart,  shown  in  Fig  II-5.  If  the  open  loop  trans¬ 
mission  of  a  unity  feedback  system  (L  *  GP,  assuming  F  »  1  for  now,  in 
Fig  II-l)  is  plotted  using  the  horizontal  and  vertical  scales  on  the 

chart,  then  at  any  given  frequency,  the  magnitude  and  phase  angle  of  T 

R 

*  L/(14L)  can  be  read  directly  from  the  curved  scales.  Conversely,  any 
point  corresponding  to  the  magnitude  and  angle  of  on  the  curved 
scales  provides  a  point  corresponding  to  the  magnitude  and  angle  of  L  on 
the  horizontal  and  vertical  scales  (Ref  3:332-334).  This  correspondence 
between  L  and  1^  on  the  Nichols  chart  is  very  important. 

Likewise,  the  Nichols  chart  can  be  used  for  the  disturbance 
response.  Recall  that  T^  ■  1/(1  +L).  By  way  of  Sidi's  transformation, 
L  ■  1/m  (Ref  1:152-155)  the  system  control  ratio  due  to  the  disturbance 
becomes  T^  *  m/(l  +m),  which  is  of  the  same  form  as  Tg  »  L/(l  +  L). 
One  could  design  the  inverse  of  the  loop  transmission,  m,  directly  on 
the  Nichols  chart,  but  it  is  much  easier  to  realize  that  by  turning  the 
Nichols  chart  upside  down,  reflecting  the  vertical  angle  of  L  lines 
about  the  -180  degree  line  (i.e.  -190  becomes  -170,  -210  becomes  -150, 
etc.),  and  reversing  the  signs  on  all  magnitude  lines,  the  chart  can  be 
used  directly  to  design  L  itself.  The  horizontal  and  vertical  lines 
still  correspond  to  the  magnitude  and  angle  of  L,  and  the  curved  magni¬ 
tude  lines  correspond  to  the  magnitude  of  (1  +  L)  (Ref  1:155).  For 
design  purposes,  only  the  magnitude  of  (1  +  L)  is  required.  Therefore 
the  curved  angle  lines  on  the  chart  can  be  ignored.  In  practice  Sidi's 


transformation  is  merely  implied  by  turning  the  Nichols  chart  upside 
down  and  modifying  the  scales  as  described  above.  The  dummy  variable. 


m,  need  not  be  considered  further. 


•.  - 

L.-  . 


Plant  Templates 

A  plant  template  is  a  plot  on  the  Nichols  chart  of  the  range  of 
uncertainty  In  the  plant  P  at  a  given  frequency  (Ref  7:290).  Consider 
the  example  P(s)  «  K/s(s+a)  where  the  gain  K  is  described  by:  2  <  K  <  8, 
and  the  location  of  the  second  pole  is  given  by:  0.5  <  a  <  2.0.  An 
infinite  number  of  possible  P's  exist  due  to  the  variation  in  parame¬ 
ters,  K  and  a;  however  each  parameter  is  a  member  of  a  set  with  finite 
boundaries.  Likewise,  the  magnitude  and  phase  angle  of  all  possible  P's 
lie  within  finite  boundaries  when  plotted  at  a  given  frequency.  The 
plant  template  is  obtained  by  plotting  Lm[P(jcj)]  vs.  Ang[P(jto)]  for  all 
possible  _P( jw)'s  at  a  given  frequency  on  the  Nichols  chart.  Note,  only 
the  outer  edges  of  the  template  need  be  calculated.  The  plant  transfer 
functions  at  the  boundaries  are  found  by  holding  one  parameter  constant 
at  a  boundary  value,  i.e.  set  K  ■  2,  and  vary  a  in  increments  from  0.5 
to  2.0.  The  frequency  response  at  to  *  1  for  the  P's  obtained  above 
provide  a  set  of  points  from  A,  (K-2,  a-0.5),  to  D,  (K-2,  a-2),  on  the 
Nichols  chart  as  shown  in  Fig  II-5.  The  process  is  continued  until  the 
complete  template  is  formed.  For  example,  for  a  -  0.5,  vary  K  from  2  to 
8  to  obtain  the  line  from  A,  (a-0.5,  K-2),  to  B,  (a-0.5,  K-8).  Tem¬ 
plates  are  needed  for  a  number  of  frequencies  taken  at  regular  inter¬ 
vals,  such  as  every  octave.  A  set  of  templates  is  shown  in  the  figure 
to  demonstrate  the  change  in  size  and  location  of  the  range  of  uncer¬ 
tainty  in  P  for  different  frequencies. 


M 
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Nominal  Plant 

To  facilitate  the  shaping  of  the  loop  transmission,  the  designer 

needs  a  reference  or  nominal  plant  transfer  function.  This  nominal 

plant,  P  ,  chosen  by  the  designer  Is  nothing  more  than  a  reference  plant 

to  be  used  in  the  definition  and  shaping  of  the  nominal  loop  transmis- 

sion,  L  ■  GP  .  There  are  no  rules  or  constraints  on  the  selection  of 
o  o 

P  .  It  doesn't  even  have  to  be  from  the  set  of  possible  P's,  but  it  is 
usually  convenient  to  choose  PQ  such  that  it  lies  at  a  recognizable 
point  on  the  templates.  It  is  convenient,  as  is  the  case  with  the 
example,  to  select  Pq  such  that  it  lies  at  the  lower,  left  corner  of  the 
templates  which  helps  keep  the  bounds  on  Lq,  to  be  described  next,  as 
near  the  center  of  the  Nichols  chart  as  possible.  Once  selected,  the  Pq 
point  should  be  marked  on  each  template,  as  in  Fig  II-5.  For  the 
example,  the  plant  described  by  PQ  -  2/(s  +  0.5)  is  chosen  as  the 
nominal  plant. 


Derivation  of  Bounds  on  L 


The  system  step  response  y(t)  is  uniquely  determined  by  the  trans¬ 
fer  function  T(s).  Likewise,  T(s),  for  a  stable,  minimum  phase  system 
(no  right-half-plane  poles  or  zeros),  is  completely  specified  by  the 
magnitude  of  the  frequency  response  |T(ju>)|  as  described  in  References  7 
and  9.  From  the  design  specifications,  the  frequency  response  of  the 
output  |Y(jw)|  can  vary  from  the  value  of  the  bound  (Ty)  to  the  value  on 
the  bound  (T^)  at  a  given  frequency  (see  Fig  II-4).  For  the  given 
example,  at  the  frequency,  1,  assume  that  |Y(jl)|  can  vary  from  0.7 
db  to  -0.8  db.  The  relative  variation  in  |Y(jl)l  is  (0.7)  db  -  (-0.8) 
db  or  1.5  db.  In  general,  the  allowable  relative  change  in  Y(jco)  at  a 


given  frequency  is  expressed  as: 


A  La  [Y(  Jo>±)  ]  -  Lm  [T^CJc^)]  -  Lm  [T^jo^)]  (II-6) 

where  T.  (Ja>)  and  T_  ( ja>)  are  the  frequency  domain  bounds  on  Y(jw). 

The  relative  change  in  the  output  is  related  to  the  control  ratio 
as  follows*  From  Fig  II-l  and  Equation  IT-1,  ALmY  ■  ALmT  •  ALm[FL/(l-+L)] 
where  L  -  GP  and  it  is  assumed  that  no  uncertainty  exists  in  G  end  F. 
Then, 


ALm  [Y( jcd^)]  -  ALm  fT( 


Ujo^) 

ALm  — - - - . 

[1  +L(jci)] 


(H-7) 


Likewise,  the  relative  change  in  L(ju>)  is  equal  to  the  relative  change 
in  the  plant. 


ALm  lUjwpj  -  ALm  [P< ja>±)]  (II-8) 

The  variation  in  P  arises  due  to  parameter  uncertainty,  thus  the  problem 
is  to  find  an  L  such  that  the  relative  change  requirements  on  the 
closed-loop  response  are  satisfied  for  the  entire  uncertainty  range  of 
P.  The  design  specifications  state  the  requirements  on  the  closed-loop 
response  Y(jcn)  and  thus  T(jui)  as  given  by  Equation  II-7.  Constraints  on 
the  loop  transmission  L(ja>)  are  desired  (Ref  7:291,  9:18). 

L(  jq>)  Bounds  on  the  Nichols  Chart 

The  relative  uncertainty  in  L  is  shown  to  be  equal  to  the  range  of 
uncertainty  in  P  by  Equation  II-8.  As  described  earlier,  the  plant 
template  is  a  plot  on  the  Nichols  chart  of  the  range  of  uncertainty  in  P 
at  a  given  frequency.  Because  Lm  (L)  ■  Lm  (P)  +  Lm  (G)  and  also 


Ang(L)  ■  Ang(P)  +Ang(G),  a  template  may  be  translated  (but  not  rotated) 
horizontally  or  vertically  on  the  Nichols  chart,  where  horizontal  and 
vertical  translations  correspond  to  the  angle  and  magnitude  requirements 
on  G(j<o)  respectively  at  a  given  frequency  (Ref  7:290).  Drawing  a  line 
on  each  of  the  templates  parallel  to  the  horizontal  or  vertical  grid 
lines  (see  Fig  II-5)  of  the  Nichols  chart  is  suggested  to  maintain 
correct  template  orientation. 

With  the  template  corresponding  to  to*  1  of  Fig  11-5,  translate  it 
to  position  1  shown  in  Fig  II-6.  Since  the  template  is  the  range  of 
uncertainty  in  P  and  L  ■  GP,  where  G  is  to  be  precisely  determined,  it 
follows  that  the  area  now  covered  by  the  template  corresponds  to  the 
variation  in  L  and  in  T  due  to  the  uncertainty  of  P.  Recall  the  corre¬ 
spondence  between  L  and  T  on  the  Nichols  chart.  Using  the  curved  magni¬ 
tude  contours,  i.e.  contours  of  constant  Lm[T(jw)],  read  the  maximum  and 
minimum  values  of  T  covered  by  the  template.  If  the  difference  between 
the  maximum  and  minimum  values  is  greater  than  the  allowable  variation 
in  T  at  the  frequency  »  1,  (  ALm  [T_(jl)J  as  given  by  Eq  II-7  and 
determined  from  Fig  II-4),  shift  the  template  vertically,  as  shown  in 
Fig  II-6,  until  the  difference  is  equal  to  ALm  [T(jl)]  (to  position 
2).  Conversely,  if  the  difference  is  less  than  that  allowed,  move  the 
template  vertically  downward  until  the  equality  is  obtained.  When  the 
position  of  the  template  achieves  the  equality  (position  2  of  the 
example),  mark  the  nominal  point  Pq  of  the  template  on  the  Nichols 
chart.  The  point  marked  corresponds  to  a  bound  on  the  magnitude  and 
phase  angle  values  of  JLq( jl )  read  from  the  horizontal  and  vertical 
scales  of  the  Nichols  chart,  where  the  nominal  loop  transmission, 
L  ( jo>. ) ,  is  given  by: 
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Repeat  the  process  horizontally  across  the  chart  at  different  values  of 
Ang(LQ).  The  points  marked  on  the  chart  form  a  curve,  repre¬ 

senting  the  boundary  of  L^jte^)  at  the  given  frequency  of  the  template. 

o^),  corre¬ 
sponding  to  w*  at  the  frequency  to*  the  variation  In  T  due  to 

the  uncertainty  In  P  is  less  than  or  equal  to  the  relative  change  in  T 
allowed  by  the  design  specifications  at  that  frequency.  Repeat  this 

o^)f  derivation  for  various  frequencies,  ^  using  the 
corresponding  plant  templates  to  obtain  a  series  of  bounds  on  L  (jcn.) 
(Ref  7:291-292). 

Likewise,  the  step  disturbance  response  specification  (line  on 
Fig  II-4),  is  converted  to  bounds  on  L^jto).  In  order  to  effectively 
reject  the  disturbance  the  following  inequality  must  be  satisfied; 

|l  +  L(jw)  I  -  l£(»l  (11-10) 

where  l£(j<o)l  is  the  magnitude  of  the  boundary,  (Tp),  on  Fig  II-4.  Con¬ 
verting  the  magnitudes  to  decibels  and  rearanging  terms,  the  inequality 
can  be  expressed  as: 

Lm  [1  +  L(jto)]  >  -  Lm  [C(jw)]  (11-11) 

Now  a  template  is  placed  on  the  inverted  Nichols  chart  such  that  its 
lowest  point  rests  directly  on  the  contour  of  constant  Lm  [1  +  L(j a>)] 
equal  to  -Lm  [C(ja>)]  at  the  frequency,  <j^,  for  which  the  template  is 
drawn.  The  point,  Pq,  is  marked  and  the  template  slid  along  *  he  same 
contour  forming  a  bound,  B_(1w. )  for  L  .  Bounds  are  formed  for  each 


boundary,  B^j 


As  long  as  L^Jfc^)  lies  outside  or  above  the  boundary,  B^Cj 


frequency,  in  this  manner  using  each  template.  Using  the  rectangular 
(Lm  L)  grid,  transcribe  the  bounds,  B^(jw^),  on  Lq  onto  the  upright 
Nichols  chart  which  already  contains  the  command  response  bounds, 
on  Lq  as  shown  in  Fig  II-7.  For  each  frequency  of  interest, 
erase  the  lower  of  the  two  Lq  bounds,  where  the  remaining  bound  is 
labeled,  The  point  here  is  that  the  worse  bound  must  be  used  in 

the  shaping  of  Lq. 

Universal  Frequency  Bound 

The  universal  frequency  (UF)  bound  ensures  the  loop  transmission, 
L,  has  positive  phase  and  gain  margins,  whose  values  depend  on  the  oval 
of  constant  magnitude  chosen  (see  Fig  II-7).  As  the  frequency,  w, 
increases,  the  plant  templates  become  narrower  and  can  be  considered 
vertical  lines  as  u  approaches  infinity.  The  allowable  variation  in  T 
Increases  with  frequency  also.  The  result  is  the  bounds  of  L^(jw^)  tend 
to  become  a  very  narrow  region  around  the  0  db,  -180  degree  point 
(origin)  of  the  Nichols  chart  at  high  frequency.  To  avoid  placing 
closed-loop  poles  near  the  j«  axis  resulting  in  oscillatory  disturbance 
response,  a  UF  bound  is  needed  on  the  Nichols  chart.  With  increasing  w, 
the  bounds  on  Lq  approximately  follow  the  ovals  encircling  the  origin. 
Choose  one  of  the  ovals  near  the  origin.  In  Fig  II-7,  the  contour  of 
constant  magnitude  equal  to  5  db  is  used  in  this  example.  From  the 
templates  corresponding  to  high  frequency,  find  the  template  with  the 
greatest  vertical  displacement,  Av,  in  db.  Av  may  be  accurately  deter¬ 
mined  by  finding  the  maximum  change  in  Lm[P(j«)j  in  the  limit  as  w  goes 
to  infinity.  Translate  the  lower  half  of  the  5  db  oval  down  the  length 
of  the  template,  i.e.  Av  db,  as  shown,  thus  obtaining  the  UF  bound  (see 


Fig  II-7)  (Note:  Professor  Horowitz  refers  to  this  bound  as  the  Univer¬ 
sal  High  Frequency  (UHF)  Bound)  (Ref  9:20-22). 


Shaping  of  the  Nominal  Loop  Transmission 

The  shaping  of  a  nominal  loop  transmission  conforming  to  the  bound¬ 
aries  of  Lq  is  a  most  crucial  step  in  the  design  process.  A  minimum 
bandwidth  design  has  the  value  of  Lq  on  its  corresponding  bound  at  each 
frequency.  In  practical  designs,  the  goal  is  to  have  the  value  of  Lq 
occurring  above  the  corresponding  bound,  but  as  close  as  possible  to 
keep  the  bandwidth  to  a  minimum.  Fig  II-7  shows  a  practical  design  for 
Lq.  Note,  any  right-half-plane  (rhp)  poles  and/or  zeros  of  Pq  must  be 
included  in  Lq  to  avoid  any  attempt  to  cancel  them  with  zeros  and/or 
poles  of  G.  Although  not  required,  using  the  poles  and  zeros  of  Pq  as  a 
starting  point  in  the  design  of  Lq  is  suggested,  to  avoid  any  implicit 
cancellation  of  roots  in  determining  G.  For  a  discussion  of  loop  trans¬ 
mission  shaping,  with  examples,  see  Appendix  A. 


Solving  for  £ 

The  compensator,  G,  is  obtained  from  the  relation:  G  -  Lo/Pq.  If 
the  Lq  found  above  does  not  contain  the  roots  of  PQ,  then  the  compen¬ 
sator  G  must  cancel  them.  Note,  cancellation  occurs  only  for  purposes 
of  design  using  the  nominal  plant  transfer  function.  In  actual  imple¬ 
mentation,  exact  cancellation  does  not  result  (nor  is  it  necessary) 
since  P  can  vary  over  the  entire  uncertainty  range. 

Provided  the  nominal  loop  transmission,  Lq,  is  shaped  properly, 
i.e.  meets  the  requirement  of  being  on  or  above  the  bound,  B^jw  ),  at 


each  corresponding  frequency,  the  variation  in  T  resulting  from  the 
uncertainty  in  P  is  guaranteed  to  be  less  than  or  equal  to  the  allowable 
relative  change  in  T  allowed  by  the  design  specifications  (Ref  7:291). 
The  design  of  the  prefilter,  F,  is  the  final  step  in  the  design  process. 

Design  of  £ 

Design  of  a  proper  Lq  only  guarantees  the  variation  in  |T(jw)|  is 
less  than  or  equal  to  that  allowed.  The  purpose  of  the  prefilter  is  to 
position  Lm[T(jw)]  within  the  frequency  domain  specifications.  For  the 
example  given  above,  the  magnitude  of  the  frequency  response  must  lie 
within  the  bounds  shown  in  Fig  11-4  which  are  redrawn  in  Fig  II-8.  One 

method  for  determining  the  bounds  on  the  prefilter,  F,  is  as  follows. 

Place  the  nominal  point  of  the  <*>■  1  template  on  the  Nichols  chart  where 
the  L^jl)  point  occurs.  Record  the  maximum  and  minimum  values  of 
Lm(T),  1.2  and  1.0  in  the  example,  obtained  from  the  curved  magnitude 
contours.  Compare  the  values  found  above  to  the  maximum  and  minimum 
values  allowed  by  the  frequency  domain  specifications  of  Fig  II-4  at  w  ■ 
1,  (0.7  db  and  -0.8  db) .  Determine  the  range,  in  db,  Lm(T)  must  be 

raised  or  lowered  to  fit  within  the  bounds  of  the  specifications.  For 
example,  at  «  *  1,  the  actual  Lm(T)  must  be  within  [Lm(Ty)  ■  0.7  db]  > 
Lm[T(jl)]  >  (Lm(TT )  -  -0.8  db] .  But,  from  the  plot  of  L  ,  the  actual 
range  of  Lm(T)  is:  1.2  db  >  Lm[T(jl)]  >  1.0  db.  To  lower  Lm[T(jl)]  from 
the  actual  range  to  the  desired  range,  the  prefilter,  Lm(F)  required  is: 

(0.7  -  1.2  db)  >  Lm[F(Jl)]  >  (-0.8  -  1.0  db),  or  -0.5  db  >  Lm[F(jl)]  >  - 

1.8  db  (See  Fig  II-8).  The  process  is  repeated  for  each  frequency  corre¬ 
sponding  to  the  templates  used  in  the  design  of  Lq.  Therefore  in  Fig 

II-9  the  difference  between  the  Ttt  and  T  curves  and  the  difference 

U  max 


between  the  and  T  ^  curves  indicate  the  requirements  for  F  as  a 

function  of  frequency. 


Fig  II-8  Requirements  on  F 

Bounds  on  F,  [LmCTy)  -  LmCT^)]  >  Lm(F)  >  [Lm(T^)  -  Lm(Tmln)]»  are 
plotted  as  a  function  of  frequency  as  shown  in  Fig  II-9.  By  use  of  the 
straight  line  approximation,  determine  a  transfer  function,  F,  such  that 
its  magnitude  lies  within  these  bounds.  The  transfer  function  obtained 
in  this  manner  is  the  prefilter,  F  (Ref  7:301). 

The  single  loop  design  is  complete  with  the  design  of  F.  The 


system  response  is  guaranteed  to  remain  within  the  bounds  of  the  design 
specifications,  provided  the  uncertainty  in  P  stays  within  the  range 
assumed  at  the  beginning  of  the  design  process  (Ref  7:288). 


III.  Multiple  Input  -  Multiple  Output  Theory 

Introduction 

The  design  approach  for  each  loop  of  the  MIMO  system  is  identical 
to  that  for  the  SISO  system  described  in  Chapter  II.  But  first  the  MIMO 
system  must  be  seperated  into  SISO  loops  which  are  equivalent  to  the 
actual  MIMO  model. 

In  general,  an  n-by-n  MIMO  system  can  be  represented  in  matrix 
notation  as  »  Pu,  where  jr  is  the  vector  of  plant  outputs,  u_  is  the 
vector  of  plant  inputs,  and  £  is  the  plant  matrix  of  transfer  functions 
relating  11  to  £.  This  £  matrix  is  formed  from  either  the  linear  differ¬ 
ential  equations  describing  the  system  or  directly  from  the  system  state 
space  representation. 

Professor  Horowitz  has  shown,  by  the  use  of  fixed  point  theory, 

that  the  Inverse  of  the  _P  matrix,  referred  to  as  (£',  contains  elements 

2 

which  are  the  inverses  of  n  single  loop  transfer  functions  equivalent 

to  the  original  MIMO  plant.  The  MIMO  problem  is  then  broken  up  into  n 

2 

loop  designs  and  n  prefilter/disturbance  problems,  which  are  each 
handled  as  described  in  Chapter  II  (Ref  10:677). 


The  MIMO  Plant 

Consider  the  multiple  input  -  multiple  output  plant  of  Fig  III-l. 
The  n  x  1  input  vector,  £,  produces  an  n  x  1  output  vector,  y_.  The 
relationship  between  j_  and  £  is  described  by  the  n  x  n  plant  matrix,  £, 
which  is  known  only  to  be  an  element  of  a  set  of  possible  £'s.  It  is 
assumed  that  the  range  of  uncertainty  in  £  can  be  determined,  probably 
in  the  form  of  empirical  data  relating  to 
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Fig  III-l  MIHO  Plant 


Note  that  the  input  and  output  vectors  are  assumed  to  be  of  the  same 
dimension.  Although  this  may  appear  to  be  a  restrictive  assumption,  it 
can  be  shown  that  with  n  inputs,  at  most  n  outputs  can  be  independently 
controlled  (Ref  6:530-536).  Thus  if  the  existing  model  defines  an 
unequal  number  of  inputs  and  outputs,  the  first  step  is  to  modify  the 
model  such  that  the  dimensions  of  the  input  and  output  vectors  are 
equal.  An  example  of  such  a  modification  is  presented  in  Chapter  V. 

The  plant  matrix,  P^,  can  be  derived  directly  from  the  set  of 
coupled,  linear,  time-invariant  differential  equations  describing  the 
behavior  of  the  plant  in  response  to  its  inputs.  Consider  a  general 
plant  model  of  the  form: 

^71  +  (bs2  +  cs)y2  -  (f)uL  +  (g)u2 

(ds)yi  +  (e)y2  -  (h)Ul  +  (i)u2  (III-l ) 

where  a  through  i  are  the  constant  coefficients,  the  y's  are  the  out¬ 
puts,  and  the  u's  are  the  inputs  to  the  plant.  The  system  of  Equation 
III-l  can  be  represented  in  matrix  notation  as; 


II-2 


/  **/ 


Although,  any  number  of  states  may  be  represented.  It  is  again 
assumed  that  the  input  and  output  vectors,  n  and  jr  respectively,  are  of 
the  same  dimension.  Assuming  the  system  is  linearized  and  the  A,  B,  and 
£  matrices  are  time  invariant,  the  plant  matrix  is; 

P  -  C[sl  -  A]-1B  (III-7 ) 

This  plant  matrix  is  actually  a  representative  member  of  a  set  of 
possible  plant  matrices  due  to  the  uncertainty  in  the  plant  parameters. 
In  practice,  a  finite  set  of  £  matrices  are  formed  representing  the 
plant  under  varying  conditions. 

MIMO  Compensation 

The  compensation  scheme  for  the  MIMO  system  is  similar  to  that  of 
the  SISO  system  of  Chapter  II.  The  basic  MIMO  control  structure  Is 
shown  In  Fig  II-3  where  £  is  the  uncertain  plant  matrix,  G  is  a  diagonal 
compensator  matrix,  and  £  is  a  prefilter  matrix.  Designs  involving  a 
non-diagonal  G  matrix  are  not  considered  in  this  thesis  (Ref  12:14). 


The  functions  of  (*  and  £  are  identical  to  those  of  G  and  F  of  the  SISO 
system  of  Chapter  II.  Fig  III-4  shows  a  more  detailed  breakdown  of  a 
two-by-two  MIMO  system  where: 


Fig  II 1-4  Two-by-Two  MIMO  System 


Constraints  on  the  Plant  Matrix 

The  set  of  J?  matrices  must  be  tested  to  ensure  that  two  critical 
conditions  are  met  (Ref  8:86-90): 

1.  £  must  not_Je  singular  for  any  possible  combination  of 
plant  parameters,  i.e.  £  must  exist. 

2.  As  s— »aB,  IPjiP22^  ^  lp12^21^  ^or  P°asible  plants. 
This  is  the  requirement  for  a  two-by  two  plant.  For  explanation  of  the 
constraint  inequality  for  the  three-by-three  or  higher  cases  see  Ref  8. 

The  first  condition  is  absolutely  necessary  to  ensure  control- 
ability  of  the  plant.  The  inverse  of  produces  the  effective  transfer 


III-5 


functions  used  in  the  design.  If  condition  2  is  not  satisfied,  it  nay  be 
possible  to  change  the  ordering  of  the  input  or  output  vector  which 
changes  the  ordering  of  the  P  matrix  elements. 


Effective  SISO  Loops 


-1  2 
Now  define  a  matrix  (J*  ■  P.  which  has  elements,  q,  The  n 


effective  transfer  functions  needed  are;  -  1/q^'. 


Ref  8  contains 


the  derivation  and  proof  of  this  equivalence.  The  n  x  n  MIMO  system  is 


now  treated  as  n  SISO  problems.  Fig  III-5  shows  the  four  effective 


SISO  loops  resulting  from  the  two-by-two  MIMO  system  (Ref  10:682). 


Each  loop  la  Fig  III-5  Is  handled  as  aa  individual  SISO  design 
problem  in  accordance  with  the  procedures  presented  in  Chapter  II.  The 
f's  and  g's  are  the  compensator  elements  F  and  G  described  previously. 


The  disturbances,  d^  represent  the  interaction  between  the  loops. 


-d  -2  — fci 

ij  k  qik 


,  k  -  i 


(III-8) 


The  b^j  in  the  above  equation  is  the  upper  response  bound,  (Ty  or 
Tp  in  Fig  II-4),  for  the  respective  input/output  relationship.  These 
are  obtained  from  the  design  specifications  (Ref  10:681-684).  Note  that 
the  first  digit  of  the  subscript  of  b^j  refers  to  the  output  and  the 
second  digit  to  the  input.  Thus  b^  is  a  function  of  the  response 
requirements  on  the  output,  y^,  due  to  the  input,  r^. 

A  recent  improvement  in  the  design  technique  involves  modification 
of  the  q's  on  the  second  and  subsequent  loops  based  on  the  g's  already 
designed.  This  reduces  the  overdesign  inherent  in  the  early  part  of  the 
design  process.  During  the  design  of  the  final  loop  the  exact  equation, 
representing  the  loop  and  the  interactions  of  the  other  loops,  is  used 
(Ref  11:977).  The  use  of  this  improvement  is  demonstrated  in  the  actual 
design.  Chapter  V. 


Basically  Non-interacting  (BNIC)  Loops 

When  the  response  of  an  output,  y^,  due  to  an  input,  r^,  is 
ideally  zero,  the  y,  .  loop  is  called  a  basically  non-interacting  (BNIC) 
loop  (Ref  10:679).  Due  to  loop  interaction  and  plant  uncertainty,  this 
ideal  response  is  not  achievable.  Therefore,  the  performance  specifica¬ 


tions  describe  maximum  responses  and  the  loop  is  handled  exclusively  as 


a  disturbance  rejection  problem 


Summary 

This  chapter  describes  the  multiple  input  -  multiple  output  plant 

*v 

and  the  plant  matrix  which  describes  it.  Guidelines  are  presented  for 
finding  the  _P  matrix,  which  relates  the  input  vector  to  the  output 
vector. 

The  division  of  the  MIMO  system  into  separate  SISO  loops  is 
presented  via  inverse  of  the  £  matrix.  After  the  equivalent  SISO  loops 
are  determined,  each  is  designed  in  accordance  with  the  SISO  design 
theory,  presented  in  Chapter  II. 


IV.  Aircraft 


Introduction 

This  thesis  concerns  a  lateral  controller  for  the  C-135  aircraft. 
A  linear  three  degree-of-freedom  lateral  model  is  generated  for  each  of 
three  flight  conditions.  The  linear,  time  invariant  (LTI)  differential 
equations  represent  the  motion  of  the  aircraft  operating  near  equilib¬ 
rium  at  each  flight  condition.  The  three  flight  conditions  chosen  are 
representative  of  most  of  the  aircraft  flight  envelope,  l.e.,  they 
represent,  as  closely  as  possible,  the  domain  of  plant  parameter  uncer¬ 
tainty.  The  changes  in  the  aircraft  flying  characteristics  from  one 
flight  condition  to  the  next  constitute  the  parameter  uncertainties 
considered  in  this  problem.  The  goal,  then,  is  to  design  a  controller 
capable  of  producing  desirable  flight  characteristics  throughout  the 
entire  flight  range  of  the  aircraft  without  the  use  of  gain  scheduling. 

The  three  flight  conditions  chosen  are  a  high  speed,  high  altitude 
cruise  of  MACH  0.75  at  42,000  feet,  a  moderate  speed,  medium  altitude 
cruise  of  MACH  0.65  at  25,000  feet,  and  a  low  speed,  low  altitude 
approach  configuration  of  275  ft/sec  at  sea  level  (Ref  5:13). 

Lateral  Equations  of  Motion 

Below  are  the  three  lateral  equations  of  motion  used  in  this 
thesis.  For  a  complete  development  of  these  equations  and  the  appli¬ 
cable  assumptions  see  Ref  5.  These  equations  are  linearized  about  a 
straight,  level,  and  unaccelerating,  equilibrium  flight  condition,  per¬ 
turbations  from  equilibrium  are  assumed  to  be  small. 


IV- 1 


Forces  in  Y  direction 


v  +  U  r  -  gttein  $  -  g$cos  0  *Yr+Yv+Y.v+Yp 
o  or  o  6  or  v  v  p 

+  Y*  8a  +  Y*  8r  (IV-1) 

IT 

Rolling  Moments 

p  -  [I  /I  lr-Lr+Lv+L,v+Lp+L{5a+L5  5r  (IV-2) 
r  1  xz  xxJ  r  v  v  pp  ®a  r 

Yawing  Moments 

r  -  [I  /I  ]p  -  N  r  +  N  v  +  N.v  +  N  p  +  8a  +  Ns  8t  (1V-3) 

*•  xz  zziV  r  v  v  pK  °a  r 

Since  perturbations  are  assumed  small,  the  following  simplifying 
assumptions  are  made.  First,  the  roll  rate,  is  assumed  to  be  equal 
to  the  angular  velocity  about  the  X-axis,  p.  Second,  the  yaw  rate,  lA, 
is  assumed  equal  to  the  angular  velocity  about  the  Z-axls,  r.  The 
actual  relationship  between  the  sideslip  angle,  0,  and  the  sideslip 
velocity,  v,  is; 


tan  0  •  v/U  (IV- 4) 

o 

The  small  angle  approximation,  tan  0  ■  sin  0  ■  0  is  assumed  valid  since 
is  small. 

These  assumptions  are  summarized  below: 


*  -  p  ( IV- 5 ) 

^  -  r  ( IV- 6 ) 


0  -  v/U 

o 


( IV- 7 ) 


Substituting  the  results  of  these  assumptions  into  the  differential 
equations  and  taking  the  Laplace  transforms  of  the  resulting  equations 


Forces  in  Y  direction 

[(1  -  Y^)s  -  Yy]0  -  [ (YpS  +  g  cos  *0)/UQ]* 

+  [(1  -  Yr/Uo)s  "  (g  8ln  *0)/Uo]* 

-  (YSa/Uo)6a  +  (Ysr/Uo)5r  (IV-8) 

Rolling  Moments 

-(L^s  +  Lq)0  +  (s2  -  Lps)# 

-  Kygs2  +  Lrs]*  -  L5a5a  +  L5r«r  (IV-9) 

Yawing  Moments 

-<y  +ya  -  +  v1* 

+  (s2  -  N  a)i  -  Nj  8a  +  Ns  8r  (IV-10) 

r  si  it 

Flight  Conditions 

The  three  flight  conditions  chosen  and  all  the  related  data  are 
obtained  from  Ref  5.  Table  IV-I  summarizes  the  flight  conditions. 


Table  IV-1 

C-135  Flight  Conditions 


Flight  Condition 

Altitude 

(ft) 

Mach  No. 

M 

Gross  Wt. 
(lbs) 

Velocity 

U  (ft/sec) 

0 

Cruise  1 

42,000 

.75 

190,000 

726 

Cruise  2 

25,000 

.65 

250,000 

660 

Power  Approach 

Sea  Level 

— 

165,000 

275 

All  angles  and  angular  rates  considered  in  this  thesis  and  imbedded 
in  the  data  are  measured  in  degrees  (Ref  3:13). 

Table  IV-2,  on  following  page,  contains  the  C-135  stability  deriva¬ 


tives  for  the  three  flight  conditions  considered. 


Table  IV-2 


C-135  Stability  Derivatives 


Parameter 

Cruise  1 

Cruise  2 

Power  Approach 

Y 

V 

-5. 7 4x1  (f 2 

-9. 46x1 O"2 

-1.2 79x1 0_1 

Y. 

V 

0 

0 

0 

Y 

P 

-1.337 

-1.583 

-2.294 

Y 

r 

2.621 

3.204 

4.277 

Y*a 

0 

0 

0 

Y*r 

16.68 

18.73 

10.55 

H 

-2.384 

-3.109 

-1.631 

H 

0 

0 

0 

L 

P 

-4.695xl0_1 

-6.381xlO_1 

-9. 074xl0_1 

Lr 

2.341X10"1 

3.248xlO_1 

5. 943x1 0-1 

* 

7.227xl(f1 

7 . 114x1 O-1 

1.433 

L*r 

2.235xl0-1 

3. 162xlO-1 

1. 223x1 0-1 

N0 

5. 089x1 0-1 

7.745xlO_1 

2.345xlO_1 

H 

l.lOxlO-2 

0 

1.6 2x1 0~2 

N 

P 

-5. 8 7x1  O'2 

-9.21xlO~2 

-1. 293x1 0_1 

N 

r 

-9. 2 7x1 0-2 

-1. 506x1 0_1 

-1. 503x1 0~l 

N*a 

3. 6 3x1 0-2 

6. OOxlO-2 

4. 03x1 O-2 

N*r 

-4.965xlO_1 

-8. 2 7 8x1 O"1 

-3. 305x1 O-1 

I 

XX 

3.602xl06 

4.013xl06 

2.813xl06 

I 

zz 

8.648xl06 

8.737xl06 

7.687xl06 

I 

xz 

-7.235xl05 

-2. 483x1 O5 

-2.561x10"* 

_V  *.• 


The  reference  axis  system  used  in  this  model  is  a  conventional 
stability  axis  system  with  the  X-axis  pointing  forward,  the  Y-axis 
pointing  to  the  right,  and  the  Z-axis  pointing  down. 

As  is  mentioned  in  the  previous  chapter,  this  model  has  three 
outputs,  0,  9,  and  and  only  two  inputs,  &a  and  *r.  In  order  to 
produce  a  square  plant  matrix,  one  of  the  outputs  must  be  eliminated. 
Since  the  sideslip  angle,  0,  and  the  yaw  angle,  4*,  are  closely  related 
(in  fact  equal  with  opposite  signs  in  the  steady  state),  the  yaw  angle 
output  is  mathematically  eliminated  from  the  model.  This  is  actually 
done  individually  for  each  flight  condition  in  the  next  chapter. 

Summary 

This  chapter  presents  the  differential  equations  and  data  that  make 
up  the  C-135  aircraft  model  used  in  this  thesis.  A  number  of  simplifi¬ 
cations  based  on  the  small  perturbation  assumption  are  Incorporated  into 
the  model. 


The  flight  conditions  chosen  represent  a  large  portion  of  the  C-135 
flight  envelope  from  a  low  altitude  approach  to  a  high  altitude  cruise. 


V.  Design 


Introduction 

This  chapter  follows  the  design  of  the  C-135  lateral  flight  con¬ 
troller  from  start  to  finish.  It  begins  with  the  manipulation  of  the 
three  degree-of-freedom  model  to  obtain  the  two-by-two  £  and  (£  matrices 
used  in  the  remainder  of  the  design  process.  Next,  the  desired  system 
output  is  modeled  to  obtain  a  frequency  domain  representation  of  the 
design  specifications.  This  is  followed  by  the  determination  of  the 
bounds  for  the  first  loop.  This  leads  to  the  shaping  of  the  first  loop 
transmission  from  which  the  first  element  of  the  compensator  matrix,  (J 
is  obtained.  Since  there  is  no  sideslip  command  to  follow,  there  are 
only  two  loops  to  consider  and  only  two  elements  of  G_  to  be  determined. 

The  recent  improvement  in  the  MIMO  design  technique,  mentioned  in 
Chapter  III,  is  used  during  the  design  of  the  second  loop.  This 
involves  using  the  exact  equation  for  the  system  during  the  design  of 
the  last  loop,  based  on  all  the  previously  designed  elements  of  (J. 

Finally,  this  chapter  presents  comments  about  the  completed  design 
and  possible  ways  to  improve  it.  In  order  to  reduce  the  large  bandwidth 
of  this  design,  a  second  design  is  presented  without  derivation.  Simu¬ 
lation  results  for  both  designs  are  included  in  Chapter  VI. 

The  Plant  Matrices 

The  first  step  in  finding  the  £  matrix  is  the  elimination  of  the 
yaw  angle  output  in  the  model  for  each  flight  condition.  Using  the 
values  of  the  parameters  for  Cruise  1  in  the  equations  for  the  model 
yields: 


V-l 


(V-1) 


(s  +  0.0574)0  +  (0.00184s  -  0.0443)*  +  (0.996s)^ 

-  (0)8a  +  (0.023)8r 

(2.384)0  +  (s2  +  0.4695s)*  +  (0.201s2  -  0.2341s)^ 

-  (0.7227)8a  +  (0.2235)8r  (V-2) 

-(0.011s  +  0.5089)0  +  (0.0837s2  +  0.0587s)* 

+  (s2  +  0.0927s )ip  -  (0. 0363)8a  -  (0.4965)8r  (V-3) 

Solving  for  ip  in  the  first  equation  gives: 

ip  -  [0.023/s]8r  -  [(1.004s  +  O.O58)/s]0 

-  [(0.002s  -  0.044)/s]*  (V-4) 

Substituting  this  into  the  other  two  equations  results  in: 

(-0.2016s2  +  0.2234s  +  2.398)0 

+  (0.9996s2  +  0.4789s  -  0.0104)* 

-  (0.7227)8a  +  (-0.0046s  +  0.2289)8r  (V-5) 

(-1.004s2  -  0.1616s  -  0.5142)0 

+  (0.0818s2  +  0.103s  +  0.0041)* 

-  (0.0363)8a  +  (-0.0231s  -  0.4986)8r  (V-6) 

These  two  equations  no  longer  explicitly  represent  the  '.raft 
rolling  and  yawing  moments  because  elements  of  the  sideforces  are  now 
included  in  each.  But  the  validity  of  the  mathematical  representation 
of  the  aircraft  in  flight  is  retained. 

To  obtain  the  matrix  representation  of  the  system,  define  an  input 
vector,  6,  and  an  output  vector,  £. 


b 


ba 

-<*>- 

t 

(V-7)  , 

I  * 

L«rJ 

-0_ 

(  V— 8 ) 


Thus  the  matrix  model  for  the  aircraft  at  this  flight  condition  is: 


0.996s2  +  0.4789s  -  0.0104  -0.2016s2  +  0.2234s  +2.398 

0.0818s2  +  0.103s  +0.0041  -1.004s2  -  0.1616s  -  0.5142. 


‘0.7727 

0.0363 


-0.0046s  +  0.2289 
-0.023s  -  0.4986 


b 


(V-9) 


For  clarity,  temporarily  name  the  two-by— two  matrices  above  M  and 

respectively.  The  goal  of  this  exercise  is  to  find  the  plant  matrix 

which  relates  the  input  vector  to  the  output  vector  by  y  *  and  the  O' 
matrix  which  is  P 

Since  My  *  N&,  then  y  *  M  yielding,  P  *  M  4  and  Q'  •  P_  ^ 

-  N-1M. 

In  general,  actual  computation  of  JP  is  recommended  for  use  in 
ensuring  the  constraints  on  it  are  met.  But,  in  the  two-by-two  case, 
the  constraints  can  be  checked  easily  using  It  should  be  obvious 

that  if  is  singular,  then  is  singular  and  doesn't  exist.  Like¬ 
wise  the  converse  is  true  if  M  is  singular.  The  first  constraint  on 

is  that  its  inverse  must  exist.  Since  neither  M  nor  N  is  singular,  both 
and  £'  exist  and  the  first  constraint  is  met. 

Using  a  computer  program,  is  computed  directly  from  M  and  _N  as 


described  above.  The  result  is: 


r 


(s  +  22.77)(s  +  0.5238) (s  -  0.0157) 
0.7278  (s  +  22.3459) 


(s  +  1 .631 )(s  ~  2.038) 
-0.051  (s  +  22.3459) 


(s  +  2.438)(s  +  0.0603) 
-0.7224  (s  +  22.3459) 


(s  +  0.087  ±  10.7945) 

0.023  (s  +  22.3459) 

(V-10) 


Notice  that  in  the  elements  of  qij%  the  constant  gain  terms 

are  displayed  in  the  denominators  and  that  the  root(s)  of  the  denomina¬ 
tors  are  identical.  As  seen  in  Chapter  III,  the  effective  SISO  transfer 
functions  used  in  the  loop  designs  are  the  inverses  of  the  elements  of 
So,  represents  the  inverse  of  the  element  in  the  ith  row  and 

jth  column  of  the  (£'  matrix  (q  ■  For  this  flight  condi¬ 

tion; 


q, ,  -  0.7278  (s  +  22.3459) 


ii 

(s  +  22.77)(s  +  0.523S)(!s  -  0.0157) 

(V— 11) 

q12  ’ 

-0.0509  (s  +  22.3459) 

(s  +  1.631)^8  -  2.038) 

(V-12) 

q21  ’ 

-0.7224  (s  +  22.3459) 

(s  +  2.438)(s  +  0.0603) 

(V-13) 

q22  " 

0.023  (s  +  22.3459) 

(s  +  0.087  ±  jO.7945) 

(V-14) 

The  second  constraint  on  jP  is  that  as  (jj  becomes  large; 
*  pllp22  *  >  *P12P21* 


(V-15) 


It's  easy  to  show  that  if  the  elements,  q^',  2."  (not  9jj)  meet  this 
constraint,  so  do  the  elements  of  P  .  In  this  case  as  co  becomes  large 
the  elements  of  in  the  above  constraint  equation  become; 


3  s  2 
s  >  s 


(V-16) 


which  is,  of  course,  true. 


The  £'s  and  the  resulting  q^  terms  for  the  other  two  flight  condi- 
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tions  are  found  in  exactly  the  same  manner.  The  results  are: 


Cruise  2 


These  two  sets  of  q's  meet  the  two  constraints  just  as  the  first 
set  did.  For  the  remainder  of  the  design  these  q's  are  used  as  the 
effective  transfer  functions  for  single  loop  designs.  Before  proceeding, 
notice  the  following  about  the  q's  for  each  flight  condition. 

The  most  obvious  thing  to  note  about  the  q's  is  that  some  of  them 
are  unstable,  or  have  poles  in  the  right  half  plane  (rhp).  This  is  not 
a  problem.  But  if  any  of  the  q's  have  zeros  in  the  rhp,  then  there  is  a 
problem.  The  solution  for  systems  with  rhp  zeros  in  the  q's  is  not 
addressed  in  this  thesis. 


Response  Modeling 

In  general,  the  design  specifications  describe  for  the  designer  the 
range  of  acceptable  responses.  To  determine  these  specifications  and/or 
to  relate  them  to  the  range  of  acceptable  responses  requires  some  engi¬ 
neering  insight  and  experience.  It  is  crucial  to  know  the  range  of  time 
responses  acceptable  for  the  system  if  Horowitz'  quantitative  method  is 
to  be  used. 

For  this  thesis,  the  design  specifications  are  left  entirely  up  to 
the  author.  The  decision  is  made  to  command  a  bank  angle  of  thirty 
degrees  and  try  to  keep  the  sideslip  angle  as  small  as  possible  through¬ 
out  the  maneuver. 

For  the  bank  angle  output,  a  first-order  step  response  model  with  a 
settling  time  of  about  four  seconds  is  considered  optimum  and  a  first- 
order  step  response  model  with  a  settling  time  of  about  ten  seconds  is 
considered  the  worst  acceptable  response.  For  the  sideslip  angle  out¬ 
put,  a  maximum  value  of  2.0  degrees  which  settles  to  zero  in  about  ten 
seconds  is  considered  the  worst  acceptable  response. 

Using  the  computer-aided-design  program,  TOTAL,  and  following  the 
procedure  described  in  Chapter  II,  transfer  functions  are  found  whose 
time  domain  responses  conform  to  the  response  descriptions  above.  The 
only  consideration  given  to  the  actual  plant  (aircraft)  model  during 
this  process  is  an  attempt  to  model  reasonable  responses. 

Appendix  B  contains  the  TOTAL  output  produced  during  the  modeling 
of  the  responses.  Included  are  printer  plots  of  the  time  and  frequency 
responses  for  each  transfer  function. 

Fig  V-l  is  the  plot  of  Lm  vs.  u>  for  the  bank  angle  response  models. 
Note  that  at  low  frequency,  there  is  very  little  difference  between  the 
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two  curves.  But  as  u>  becomes  larger,  the  difference  between  them 
Increases.  Had  the  specifications  allowed  for  some  overshoot  in  the 
optimum  response  the  difference  between  the  two  curves  would  become 
larger  at  low  frequency. 


Fig  V-l  Bounds  on  Bank  Angle  Response 


In  this  figure,  the  b^  and  a^  curves  correspond  directly  to  the  T  and 
curves  of  Fig  II-4  respectively.  The  subscripts  indicate  that  they 
are  bounds  on  the  first  output  response  due  to  the  first  command  input. 

Fig  V-2  is  the  plot  of  Lm  vs.  O)  for  the  sideslip  angle  response 
model.  Note  that  its  magnitude  never  gets  above  zero.  This  is  typical 
of  a  basically  noninteracting  (BNIC)  loop  in  a  MIMO  control  system.  The 
b£^  of  Fig  V-2  corresponds  to  the  TQ  curve  of  Fig  II-4.  The  subscript 


indicates  that  this  is  the  bound  on  the  second  output  response  due  to 
the  first  command  input. 


Fig  V-2  Bound  on  Sideslip  Angle  Response 

These  plots  of  response  model  magnitudes  vs.  frequency,  based  on 
the  design  specifications  are  all  the  information  needed  to  Insure  the 
final  system  responses  are  acceptable.  These  plots  are  used  to  form  the 
boundaries  for  the  transfer  functions  of  the  effective  SISO  control 
loops  during  the  design. 

Selection  of  the  First  Loop  to  Design 

In  this  problem  there  is  no  sideslip  command,  therefore,  only 

two  of  the  four  SISO  loops  of  Fig  III-5,  resulting  from  the  two-by-two 
MIMO  system  need  to  be  considered.  These  are  the  bank  angle  output  due 
to  a  bank  command,  y^,  and  the  sideslip  angle  output  due  to  a  bank 
Fig  V-3  is  Fig  III-5  redrawn.  Recall  that  d^  and  d^2 


command,  y91. 


Fig  V-3  Effective  SISO  Loops 


Of  the  two  remaining  loops,  one  must  be  selected  for  consideration 
first.  There  is  no  hard  and  fast  rule  as  to  which  loop  to  design  first 
but  in  general  BNIC  loops  are  easier  since  there  is  no  command  to 
consider.  Notice  that  in  the  1°°P  of  Fig  V-3  the  effects  of  input, 
r^,  can  follow  either  of  two  paths  to  the  output,  y^l*  enters  as 
part  of  the  disturbance,  ^21*  an<*  directly  through  fj^.  Since  it  is 
deslreable  for  ^  t0  ^e  zero,  the  less  effect  r^  has  on  the  loop,  the 
better.  Nothing  can  be  done  to  eliminate  the  disturbance,  but  f£^  is  to 
be  designed  and  can  be  made  zero,  leaving  only  the  disturbance  to  con¬ 
sider.  The  second  loop,  the  sideslip  output  due  to  a  bank  command  is  a 
BNIC  loop  and  is,  therefore,  designed  first.  Thus,  one  is  left  with  a 
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simple  disturbance  rejection  problem.  Fig  V-4  shows  the  loop  as  it 
is  dealt  with  for  this  problem. 


Fig  V-4  Diagram  of  y 21  Loop 


Bounds  on  BNIC  Loop  Transmission 

Recall  that  the  disturbance  term,  d^^,  is  I ^ I / I 1 1  where  the  b^ 
term  is  the  upper  of  the  two  bank  response  curves.  This  provides  a 
worst  case  approximation  for  the  disturbance. 

In  order  to  ensure  rejection  of  the  disturbance,  the  following 
inequality  must  be  satisfied: 


the  above  inequality.  At  any  given  frequency,  the  right  side  of  the 
inequality  reduces  simply  to  a  number  greater  than  zero.  It  is  helpful 
to  have  Bode  plots  of  each  of  the  terms  on  the  right  side.  With  these, 
one  can  simply  add  and  subtract  the  magnitudes  at  each  frequency  to 
obtain  a  minimum  for  Lm  (1  +  12)  in  db.  For  clarity,  define  a  new 
variable,  E(u>),  such  that; 


thus; 


E(U» 


b1l(W)  l£22(ju»)l 

b21(U,) 


(V-27) 


Lm  (E)  -  Lm  (bu)  +  Lm  (q22>  -  Lm  (b21)  -  Lm  (qn)  (V-28) 

By  way  of  the  transformation,  m  -  1/12,  turn  the  Nichols  chart 

upside  down,  as  described  in  Chapter  II.  Lm  (1  +  1^  >  -3  db  is  chosen 

as  the  UF  bound,  based  on  phase  margin  considerations  [remember  to 

reverse  the  signs  on  the  Lm  (1  +  12)  curves]  (see  Eqs  Il-iO  and  II).  At 

this  point,  there  are  a  number  of  ways  to  obtain  the  bounds  on  1.  .  If 

IQ 

^22o  *3  c*losen  8uch  that  It  lies  at  the  lowest  point  of  the  templates  at 
all  frequencies,  then  the  bounds  on  Lm  (1  +  12)  which  are  greater  than  0 
db  can  be  used  directly.  Simply  select  the  flight  condition  which 
produces  the  largest  Lm  [E(w^)]  and  draw  in  the  it)  ■  bound  along  the 
curved  Lm  (1  +  1^)  line  corresponding  to  the  value  of  Lm  [E(u>^)], 

realizing  that  any  possible  q22  other  that  q22(j  puts  Lm  (12)  and  the 
resulting  Lm  (1  +  12)  above  the  boundary  at 

When  the  minimum  value  for  Lm  (1  +  12)  is  less  than  the  UF  bound, 
the  process  becomes  a  little  more  complicated.  And,  since  the  magni¬ 
tudes  for  the  three  q22's  vs*  frequency  for  this  problem  cross  each 


other  (i.e.  no  lowest  ^220  can  chosen),  another  scheme  must  be  found. 
Consider  the  following: 


12o  "  g2q22o  * 


82q22q22o 

q22 


12q22o 

q22 


(V-29) 


Using  the  Nichols  chart,  the  minimum  for  Lm  (1  +  12)  at  each  flight 
condition  is  converted  into  a  minimum  for  Lm  (1^)  at  regular  phase  angle 
intervals.  When  the  Lm  [E(u^)j  is  less  than  -3  db  (UF  bound),  then  -3 
db  becomes  the  bound.  Since  the  values  are  in  db,  the  multiplication 
and  division  indicated  in  Equation  V-29  become  addition  and  subtraction. 
When  the  magnitude  and  angle  differences  between  ^220  an<*  q22  at  a 
particular  flight  condition  and  frequency  are  added  to  the  minimum  for 
1^,  the  result  is  a  minimum  for  the  magnitude  and  angle  of  1 at  that 
flight  condition  and  frequency.  This  is  done  for  each  of  the  three 
flight  conditions  and  the  worst  case  boundary  of  the  three  is  chosen  as 
the  boundary  for  l2Q-  Below  are  two  examples  of  this  process. 

Example  1_  The  9220  chosen  for  this  problem  is; 


q22o 


0.03(8  +  30) 

(s  +  0.15  +  jl) 


(V-30) 


Notice  that  this  is  not  one  of  the  actual  922~s  as  liste<*  *n  Equations 
(V-14),  (V-20),  or  (V-24).  There  is  no  requirement  for  ^220  t0  an 
actual  q  of  the  system.  Because  no  convenient  ^220  ex*sts  ^ or  this 
problem,  a  simplified  form  of  for  the  flight  condition.  Cruise  2,  is 
chosen  as  922o' 

Table  V  -  1  contains  the  relevent  information  for  the  point,  where 
0.1  rad/sec.  The  Information  is  obtained  from  the  Bode  plots  of  the 
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responses  and  of  the  q's. 


Table  V-l 

Cruise  1 

Cruise  2 

Pow.  App. 

,blll 

- > 

-0.17  db 

< - 

1  b2 1 1 

- > 

-52.6  db 

< - 

q22o 

- > 

-1.0  db  -1.5°  < - 

q21 

35.1  db 
-241° 

36.6  db 
-223° 

14.5  db 
-229° 

q22 

-1.8  db 
-1.3° 

-1.7  db 
-1.3° 

-6.1  db 
-4.1° 

In  order  to  satisfy  Equation  V-26,  Lm  (1  +  for  all  three  flight 
conditions  at  U)  m  0.1  rad/sec  must  be  greater  than  Lm  (E).  This  works 
out  to  15.5  db  for  Cruise  1,  14.1  db  for  Cruise  2,  and  31.8  db  for  the 
Powered  Approach.  The  Powered  Approach  condition  produces  the  largest 
bound  and  therefore  dominates  at  this  frequency.  On  the  inverted 
Nichols  chart,  note  that  the  Lm  (1  +  1)  ■  32  db  curve  is  nearly  a 
horizontal  line,  which  deviates  from  the  Lm  (1)  -  32  db  line  very 
little.  This  makes  sense,  since  1  »  1  ,  Lm  (1  +  1)  ■  Lm  (1)  is  a  good 
approximation.  The  idea  here  is  not  to  make  the  approximation,  but  to 
point  out  the  relationship  between  Lm  (1)  and  Lm  (1  +1)  on  the  Nichols 
chart  more  clearly. 

For  the  Powered  Approach  condition,  the  difference  between  q£20  an<* 
q22  is  db  and  2.6°.  Adding  these  values  to  the  bound  of  Lm  (1  +  1 ,, ) 
■  31.8  db  results  in  a  curve  that  has  the  same  shape  as  Lm  (1  +  1)  - 
31.8  but  shifted  5.1  db  up  and  2.6°  to  the  right.  This  curve  is  now  the 
bound  on  !  1  +  lj  I  and  on  for  this  flight  condition. 


This  is  repeated  for  the  other  two  flight  conditions  but  at  0.1 


rad/sec  the  other  two  flight  conditions  produce  boundaries  that  are 
below  the  Powered  Approach  bound.  Therefore  the  Powered  Approach  bound 
dominates  at  W  -  0.1  rad/sec  and  is  the  one  used  during  the  design  of 


Example  2  At  U)m  2  rad/sec  the  process  is  repeated  to  demonstrate 
how  to  deal  with  the  UF  boundary.  Table  V-2  contains  the  relevent 
information  for  W*  2  rad/sec. 


Table  V-2 


Cruise  1 

Cruise  2 

Pow.  App. 

Ih  1 

-12.5  db 

bll 

v.”- 

Ih  1 

. 

-25.6  db 

"21 

.tn  ;  jk  _u;° 

q22o 

J  QD 

q21 

8.2  db 
-303° 

9.2  db 
-212° 

-7.9  db 
-240° 

q22 

-16.3  db 
-169° 

-10.8  db 
-166° 

-26.5  db 
-145° 

Again  the  value  of  Lm  (E)  is  the  maximum  for  Lm  (1  +  l^).  This  is 
-11.4  db  for  Cruise  1,  -6.9  db  for  Cruise  2,  and  -5.5  db  for  the  Powered 
Approach  condition.  All  three  of  these  values  are  less  than  -3  db. 
Therefore  the  UF  boundary  value  of  -3  db  prevails  at  this  frequency  and 
(as  it  works  out)  all  frequencies  higher  than  2  rad/sec. 

So,  for  all  three  flight  conditions  the  boundary  for  Lm  (1  +1^)  is 
the  oval  corresponding  to  -3  db  on  the  inverted  Nichols  chart.  To  get 
the  bounds  on  Lm  (1  +  ^  ),  rePeat  the  steps  taken  for  UJ-  0.1  rad/sec. 
For  example,  for  Cruise  1,  the  difference  between  ^220  and  q22  is  db 
and  4°.  The  bound  for  Cruise  1  is  the  exact  shape  of  the  -3  db  oval 
shifted  5.8  db  up  and  4°  to  the  right.  If  this  is  repeated  for  the 


ocher  two  flighc  conditions,  each  shifts  the  oval  a  different  amount. 
The  overall  boundary  on  |l  +  1.  i  and  thus  on  1-  is  the  outer  boundary 

40  40 

of  the  shape  formed  by  this  shifting  of  the  -3  db  oval  for  the  three 
flight  contidions. 

Bounds,  Bq(Ju^),  are  found  in  this  manner  at  frequency  intervals  of 
about  every  octave  between  w*  0.1  and  u)*  500  rad/sec.  Fig  V-5  shows 
the  bounds  for  1.  . 

40 


Shaping  of  BNIC  Loop  Transmission 

Once  the  bounds  on  l2o  are  determined,  the  loop  transmission  is 
shaped  as  described  in  Chapter  II.  In  this  case  the  dominate  bound  is 
at  u)  ■  1  rad/sec.  In  other  words,  if  Lm  (12q)  is  at  or  greater  than 
about  36  db  at  U  ■  I,  all  the  bounds  at  lower  frequencies  are  also  met 
(assuming  a  relatively  constant  slope  of  -6  db/decade  due  to  a  pole  at  u) 
»  0  rad/sec).  Therefore,  Lm(l2o>  equals  36  db  at  u*  1  rad/sec  is  used 
as  a  starting  point  in  the  shaping  of  12q. 

At  and  beyond  W ■  2  rad/sec  the  only  consideration  is  ensuring  the 
loop  transmission  does  not  enter  the  UF  bounds  which  are  represented  by 
the  ovate  shapes  on  the  Nichols  chart.  Fig  V-6  shows  the  Nichols  chart 
with  the  boundaries,  B^(juJ!),  and  the  shaped  12q. 

The  final  mathematical  description  of  the  12q  chosen  is: 


58.2  (l  +  s)2(l  +4>(1  +  |^) 


i(l  +2)(1  +2o)fl  +  (1.01)2  +(l.0lJJfl  +  (190)2  +(l9(J1 

( V— 3 1 ) 


Dividing  through  by  q__  to  get  g7  yields; 


66.1  (1  +  s)  (1  +  ■=■) 


s(l  +  ■>)(!  +  +  (190)^  ^ 


This  loop  is  driven  entirely  by  the  effective  disturbance, 
fore,  there  is  no  command  and  no  F.  matrix  element  for  this  loop 
this  element  of  £  designed  it's  time  to  turn  to  the  other  loop. 


Effective  Plant  for  Second  Loo 


The  entire  compensated  system  can  be  represented  by  a  signal  flow 
diagram  as  shown  in  Fig  V-7.  The  p^'s  are  the  elements  of  the  plant 
matrix,  P. 


Fig  V-7  Compensated  System 


Note  that  Fig  V-7  is  a  simplification  of  Fig  III-4.  By  use  of  Mason's 
Gain  Formula  (Ref  3:162-164)  the  transfer  function,  t^  ■  y^/r^,  can  be 
obtained.  The  following  substitutions  are  used  to  obtain  t^  in  terms  of 
the  q^'s  instead  of  the  p^j's: 


det(P) 


P11P22  "  P12P21 


(V-33) 
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From  Fig  V-7  t^  is: 

fll  81  P11  (1  +  g2  P22^  “  fl  81  P21  g2  p12 

C11 

1  +  gl  Pu  +  g2  P22  “  81  P21  g2  Pi2  +  gl  P11  g2  P22  (V-40) 

Rearanging  terms  and  substituting  for  the  p^'s  yields: 

fll  gl  t  CA/q22)  +  82 A  ] 

cn  *  - 

1  +  (gL  A)/q 22  +  (g2  A)/qu  +  g:  g2A  (V— 41) 

Dividing  through  by  A  and  substituting  Equation  V-39  for  1/A  yields: 

_ fll  gl  ^q22  +  g2^ _ 

a  1/(qll  q22J  "  1/(q12  q21)  +  gl/q22  +g2/qll  +  gl  g2  (v_42) 

Now  multiply  through  by  q^q22  and  substitute  Equation  V-38  into  the 
denominator  to  get: 


t  _  fU  fl  (qU  +*2  qll  q22> _ 

11  1  -  Y  *  *1  In  +  *1  «u  *2  q22  +  *2  q22  <v'43> 

Combining  terms  and  using  Equation  V-37,  the  resulting  actual  equation 
for  t^  is; 


f-LlV1  +12> 


(1  +  1^(1  +  12)  -  Y 


(V-44) 


Now  if  a  new  loop  transmission,  1^  is  defined  as; 


1  +  In 


(V-45) 


then  t^ ^  can  be  manipulated  into  (Ref  12:102); 


fll  X1 


1  +  1, 
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Now  1, 


*  * 
g.q..  where  q. .  is; 


1  +  1, 


(V-47) 
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Since  g2  has  already  been  designed,  12  is  known  and  can  be  deter¬ 


mined.  Using  TOTAL,  q^  is  found  for  each  flight  condition.  The 


resulting  transfer  functions  can  be  found  in  Appendix  C.  Fig  V-8  shows 


the  Bode  plots  of  q^  for  all  three  flight  conditions. 


The  design  of  1^  can  now  be  accomplished  as  a  simple  single  loop 


*  * 
problem  where  q^  is  the  nominal  plant  transfer  function.  The  q^ 


for  flight  condition,  Cruise  2,  is  chosen  as  q^  since  it  has  the 


smallest  magnitude  for  all  frequencies.  Bounds  on  1^  are  found  based 


on  the  allowable  variation  in  t..  from  the  response  plots,  Fig  V-l.  Fig 


V-9  is  the  Nichols  chart  with  templates  of  q^  and  bounds  on  1^q 
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Fig  V-8  Plots  of 


Shaping  of  the  Bank  Response  Loop 

a 

With  the  bounds  on  1^q  defined,  the  loop  transmission  is  shaped 
exactly  as  described  in  Chapter  II.  Although  it  is  not  as  obvious  as  in 
the  last  loop,  the  dominant  bound  occurs  at  u»  0.4  rad/sec.  Fig  V-10 

A 

is  the  Nichols  plot  of  the  resulting  1^q  which  is: 


-43.4  (1  +  q~)(1  +  s) 
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A  A 

Notice  the  unstable  pole  of  1^q  .  This  pole  is  part  of  q^Q  and 


it  is  very  important  that  any  unstable  poles  of  the  nominal  q 


11 


be 


included  in  lj.io  *  If  the  loop  is  shaped  without  including  an  unstable 
pole,  solving  for  g^  results  in  a  rhp  zero  cancelling  it.  Since  one 
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cannot  rely  on  exact  cancellation,  then  an  unstable  system  mey  result. 
If  the  rhp  poles  are  included,  the  loop  can  be  shaped  to  ensure  adequate 


stability  and  phase  margins,  as  is  the  case  here.  The  resulting  g^  is; 


1.29  (  1  +  s) 

'l +T5T4 +(f)2'2  (v'49) 


Prefilter  Design 

There  is  only  one  prefilter  element  to  be  designed  for  this  prob¬ 
lem.  That  element  is  f^.  As  explained  in  Chapter  II,  ensures  that 
the  variation  in  the  magnitude  of  t^,  as  a  function  of  frequency, 
remains  within  limits  specified  by  the  difference  between  the  magnitudes 
of  the  allowable  responses.  The  prefilter  element  f^  is  to  ensure  that 
the  actual  magnitude  of  t^  remains  between  those  bounds. 

Fig  V— 11  is  the  Bode  plot  of  the  three  t^'s  resulting  from  the 
three  flight  conditions  superimposed  on  the  t^  response  bounds.  This 
plot  indicates  the  form  the  frequency  response  of  f^  must  take. 

The  prefilter  element  chosen  is  of  the  simplest  form  possible. 
Note  that  this  f^  only  ensures  the  proper  t^  in  the  low  frequency 
region.  A  more  complicated  f^  may  ensure  a  proper  tj^  for  a  larger 
frequency  range  but  a  more  complicated  f^  is  not  necessary  to  meet  the 
specifications.  The  f^  chosen  is; 

f,,  "  —0.7  -  1 

(1  +  Q^7)  (V-50) 
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Fig  V— 1  1  Bounds  and  Resulting  t^'s 

This  completes  the  design  of  the  compensator  and,  as  shown  in  the 
next  chapter,  it  meets  the  specifications  easily.  The  bandwidth  of  the 
g2  element  is  quite  large  and  may  be  unreasonable  in  terms  of  suscepti¬ 
bility  to  sensor  noise  and/or  implementation  of  the  element  transfer 
function. 


Alternative  Design 


In 

an  effort 

to  reduce 

the  bandwidth  of  g2  a 

second 

design 

is 

offered. 

Remember 

that  the 

dominate  bound  on 

occurs 

at  W  » 

1 

rad/sec  and  is  nearly  constant  at  36  db.  The  difference  between  the 
straight  line  approximation  of  the  magnitude  of  and  the  exact  value 
at  <j)  m  1  rad/sec  is  about  15  db.  This  is  due  to  the  very  small  damping 


factor  of  the  <*22o  P°*e  Pa*r  at  1*  Notice  that  each  ^as  a 
similar,  underdamped  pole  pair.  In  order  to  take  advantage  of  this, 
another  12q  is  shaped  ignoring  the  bound  at  u>»  1. 

The  compensator  element  g2  f°r  c^e  second  design  is; 


14-4<1+7072  +(o)  )(1+f> 

,a  +sKi  +fj)ti  +^2  +  (fr)2]2 


(V-51) 


Since  g2  *-s  different,  the  effective  plant,  q^  is  different  also. 

Therefore  the  entire  design  process  for  g^  is  reaccomplished  beginning 

* 

with  the  computation  of  for  each  flight  condition.  In  this  case, 

* 

by  coincidence,  the  frequency  characteristics  of  the  new  q^  ' s  change 
very  little  as  a  result  of  changing  g2«  The  gL  and  the  fn  of  the  first 
design  also  meet  the  respective  requirements  for  the  second  design  and 
are  left  unchanged. 

It  will  be  seen  in  Chapter  VI  that  the  results  for  the  second 
design  are  not  acceptable.  It  is  presented  here  for  comparison  and  to 
demonstrate  the  consequences  of  ignoring  a  boundary. 


Summary 


This  chapter  follows  the  design  of  the  three  required  controller 
elements,  g2,  g^,  and  f^,  from  the  elimination  of  the  yaw  angle  output 
from  the  model,  through  the  actual  formulation  of  the  compensator  trans¬ 
fer  functions.  Because  of  the  large  bandwidth  of  the  g2  element  of  the 
first  design,  a  second  design  is  also  presented. 

The  results  of  computer  simulations  for  both  designs  are  presented 


in  Chapter  VI.  A  discussion  of  these  results  is  presented  in  Chapter 
VII. 
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VI 


Simulation  Results 


Introduction 

This  chapter  presents  the  results  of  simulations  for  both  designs 
presented  in  Chapter  V.  All  the  simulations  are  products  of  the  com¬ 
puter-aided-design  program,  TOTAL.  The  input  for  each  simulation  is  a 
step  command  input  demanding  thirty  degrees  of  bank. 

The  first  design  meets  the  specifications  set  forth  in  Chapter  V 
easily.  The  second,  which  is  the  result  of  ingnoring  the  dominant  bound 
on  the  second  loop,  does  not  meet  the  specifications.  It  is  presented 
here  as  an  example  of  the  consequences  of  ignoring  a  boundary.  Transfer 
functions  for  all  simulations  are  in  Appendix  C. 

Design  One 

Figs  VI-1  through  VI-4  are  the  results  of  simulations  for  the  high 
altitude,  high  airspeed  condition.  Cruise  1.  Fig  VI-1  is  the  bank  angle 
response  to  a  thirty,  degree  step  bank  angle  command.  The  desired  first- 
order  response  is  achieved  with  about  a  four  second  settling  time.  Fig 
VI-2  is  the  deflection  of  the  control  wheel  during  the  maneuver.  Fig 
VI-3  Is  the  sideslip  response.  Note  that  the  peak  value  of  sideslip  is 
less  than  0.02  degrees  throughout  the  maneuver.  Fig  VI-4  is  the  deflec¬ 
tion  of  the  rudder  during  the  maneuver.  As  expected,  the  aircraft 
response  at  this  flight  condition  is  excellent. 


VI-1 
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Fig  VI-1  Bank  Angle  Response  -  Cruise  1  -  Design  1 


Fig  VI-2  Control  Wheel  Deflection 


Cruise  1 


Design  1 


Fig  VI-4  Rudder  Deflection 


Cruise  1 


Design  1 


Figs  VI-5  through  VI-8  are  the  simulation  results  for  the  medium 
altitude,  moderate  airspeed  flight  condition,  Cruise  2,  to  the  same 
thirty  degree  of  bank  command  input.  Although  the  controller  is 


unchanged,  the  responses  are,  again,  acceptable. 

Fig  VI-5  is  the  bank  angle  response.  The  desired  first-order 
response  is  again  achieved  but  the  settling  time  is  about  six  seconds. 
Fig  VI-6  is  the  control  wheel  deflection  during  the  maneuver.  Its 
response  is  very  similar  to  that  of  Cruise  1.  Fig  VI-7  is  the  sideslip 
angle  response.  The  peak  value  is  on  the  order  of  0.01  degrees  which  is 
even  less  than  that  of  Cruise  1.  Fig  VI-8  is  the  rudder  deflection 
during  the  maneuver. 


RANK  ANGLE  RESPONSE 


CRUISE  2 


Fig  VI-8  Rudder  Deflection  -  Cruise  2  -  Design  1 

Figs  VI-9  through  VI-12  are  the  aircraft  responses  for  the  low 
altitude,  low  airspeed  flight  condition  in  the  approach  configuration, 
Powered  Approach.  As  expected,  the  aircraft  response  is  degraded  some¬ 
what  but  the  performance  is  still  well  within  the  specifications  of 
Chapter  V. 

Fig  VI-9  is  the  bank  angle  response.  It  is  still  first-order  but 
the  settling  time  is  about  seven  seconds.  Fig  VI-10  is  the  control 
wheel  deflection  during  the  maneuver.  Fig  VI-11  is  the  sideslip  angle 
response.  The  peak  value  is  about  0.15  degrees,  which  is  about  a  factor 
of  ten  worse  than  that  of  Cruise  1,  but  still  well  within  specifica¬ 
tions.  Fig  VI-12  is  the  rudder  deflection  during  the  maneuver.  Notice 
that  the  response  is  much  different  than  that  of  either  Cruise  1  or  2. 
This  is  due  to  the  change  in  spoiler  operation  in  this  configuration. 


Fig  VI-9  Bank  Angle  Response 
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Fig  VI-10  Control  Wheel  Deflection 
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The  three  flight  conditions  are  chosen  to  represent  a  major  portion 
of  the  0135  flight  envelope  and  are  limited  by  the  availability  of 


flight  data.  These  simulation  show  that  the  controller  produces  accept¬ 
able  aircraft  responses  at  all  three  flight  conditions.  It  follows 
logically  that  the  aircraft  performance  at  any  flight  condition  between 
these  extremes  would  also  be  acceptable. 

Design  Two 

Figs  VI-13  and  VI-14  are  the  results  of  simulations  for  the  high 
altitude,  high  airspeed  condition,  Cruise  1.  Fig  VI-13  is  the  bank 
angle  response  to  a  thirty  degree  bank  angle  step  command.  This 
response  is  virtually  unchanged  from  that  of  Design  One.  Fig  VI-14  is 
the  sideslip  angle  response  to  the  same  command.  The  maximum  value 
specification  is  met  but  the  settling  time  specification  is  not. 


Fig  VI-14  Sideslip  Angle  Response  -  Cruise  1  -  Design  2 

Figs  VI-15  and  Vl-16  are  the  simulation  results  for  the  medium 
altitude,  moderate  airspeed  flight  condition.  Cruise  2.  Fig  VI-15  is 
the  bank  angle  response,  which  is  nearly  identical  to  the  response  of 
Design  1  and  still  well  within  design  specifications.  Fig  VI-16  is  the 
sideslip  angle  response.  Again,  the  peak  value  specification  is  met, 
but  the  settling  time  is  far  too  long. 


Fig  VI-15  Bank  Angle  Response  -  Cruise  2  -  Design  2 
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Fig  VI-16  Sideslip  Angle  Response 
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Summary 

This  chapter  presents  the  results  of  simulations  for  both  the 
controllers  designed  in  this  thesis.  The  first  controller,  designed  in 
full  compliance  with  Dr.  Horowitz'  technique,  meets  the  specifications 
over  the  entire  range  of  uncertainty  considered.  The  second  design, 
which  is  the  result  of  an  intentional  violation  of  a  dominant  boundary, 
produces  results  which  are  within  the  peak  value  specifications  but 


which  violate  the  settling  time  requirements. 


VII.  Conclusions  and  Recommendations 


Thesis  Summary 

This  thesis  is  intended  to  describe  and  demonstrate  the  use  of  Dr. 
Horowitz'  Quantitative  Feedback  technique  to  design  a  lateral  controller 
for  a  large  transport  aircraft,  specifically  the  C-135.  The  first 
controller  presented  easily  meets  all  the  design  specifications  over  the 
entire  range  of  uncertainty  considered.  The  fact  that  the  maximum  peak 
sideslip  angle  is  more  than  ten  times  less  than  that  allowed  by  the 
specifications  is  significant  and  is  dealt  with  in  the  discussion  sec¬ 
tion  of  this  chapter.  In  addition  to  conclusions  and  recommendations, 
this  chapter  presents  a  brief  comparison  of  this  design  to  that  of  a 
concurrently  produced  thesis  concerning  the  use  of  Dr.  Brian  Porter's 
method  to  design  a  controller  for  the  C-135  aircraft. 

Chapter  I  presents  an  introduction  to  the  C-135  lateral  flight 
control  problem.  Chapters  II  and  III  present  an  overview  of  Quantita¬ 
tive  Feedbeck  Synthesis  theory.  The  single  input  -  single  output  prob¬ 
lem  is  described  in  Chapter  II  and  the  extension  to  the  multiple  input  - 
—multiple  output  problem  is  presented  in  Chapter  III.  Chapter  IV 
describes  the  aircraft  model  used  in  the  design.  The  actual  design  of 
the  controller  is  presented  in  Chapter  V.  The  simulation  results  for 
both  the  primary  design  and  an  alternate  design  are  shown  in  Chapter  VI. 

Appendices  to  this  thesis  contain  examples  of  loop  transmission 
shaping  (Appendix  A)  which  are  provided  courtesy  of  Dr.  Horowitz, 
response  models  generated  during  the  design  (Appendix  B) ,  and  the  actual 
transfer  functions  of  the  effective  plants  and  of  the  compensated  system 
(Appendix  C). 


Discussion 


The  specifications  on  which  both  designs  are  based  allow  a  peak 
value  of  2.0  degrees  of  sideslip  in  response  to  a  30  degree  bank  com¬ 
mand.  The  worst  sideslip  response,  produced  by  the  Power  Approach 
flight  condition,  reaches  a  peak  value  of  only  about  0.15  degrees.  This 
is  only  about  7.5%  of  the  allowable  value.  One's  initial  reaction  to 
learning  that  the  design  "beats"  the  specifications  by  so  much  might  be 
delight,  but  it  is  an  indication  of  serious  overdesign  and  a  waste  of 
compensation.  Although  the  problem  is  not  solved  here,  this  section 
deals  with  the  possible  causes  and  provides  a  suggested  approach  to 
reduce  the  overdesign. 

The  first  and  most  obvious  place  to  look  for  the  overdesign  is  in 
the  sideslip  loop.  It  is  this  loop  which  controls  the  sideslip  angle 
and  as  noted  previously  the  compensator  in  the  sideslip  loop  has  a 
rather  large  bandwidth.  Inspection  of  the  bound  derivation  and  actual 
loop  design  reveals  little  overdesign.  The  combination  of  low  frequency 
bounds  and  stability  considerations  dictate  a  loop  transmission  with  a 
bandwidth  very  close  to  that  of  Design  1.  Some  bandwidth  (and  therefore 
some  compensation)  reduction  could  be  achieved  but  not  enough  to  account 
for  the  magnitude  of  the  overcompensation. 

A  comparison  between  the  response  model  used  in  the  design  of  the 
sideslip  loop  and  the  sideslip  responses  from  the  simulations  shows  that 
the  aircraft  reacts  much  more  slowly  than  assumed  during  the  modeling 
process.  The  resulting  design  appears  to  provide  compensation  at  much 
higher  frequency  than  necessary.  In  other  words,  the  response  far 
exceeds  the  design  specifications  because  the  response  model  demanded  a 
compensator  with  a  much  larger  bandwidth  than  needed  to  control  the 
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aircraft. 


Although,  time  does  not  permit  a  thourough  investigation  here,  the 
bandwidth  of  the  sideslip  compensator  could  be  reduced  greatly  by  re¬ 
shaping  the  sideslip  loop  transmission  based  on  a  slower  response  model. 

Comparison  to  Porter  Design 

Capt  William  Locken,  GE-83D,  has  written  a  thesis  concerning  flight 
control  laws  for  the  C-135  using  the  method  of  Professor  Brian  Porter 
(Ref  13).  This  section  is  intended  as  a  qualitative  comparison  between 
the  Porter  technique  and  the  Horowitz  technique.  By  necessity,  much  of 
this  section  is  personal  opinion  on  the  part  of  the  author. 

The  primary  difference  between  the  two  methods  is  that  the  Horowitz 
approach  treats  plant  uncertainty  (robustness)  and  performance  specifi¬ 
cations  as  given  information  on  which  the  design  is  based,  while  the 
Porter  approach  involves  testing  for  robustness  and  conformance  with 
specifications  after  a  proposed  design  is  completed. 

A  Porter  design  can  be  completed  much  more  quickly  than  a  Horowitz 
design,  if  performance  tolerances  are  met  in  the  first  few  iterations. 
It  is  possible,  however,  that  a  Porter  design  could  take  considerably 
longer  to  complete  if  the  specifications  are  difficult  or  impossible  to 
meet.  The  Horowitz  approach  provides  insight  into  the  design  problem 
and  difficult  or  impossible  specifications  can  be  easily  identified, 
giving  the  de'-'gner  the  information  needed  to  modify  r  e  tolerances  or 
abandon  the  problem  if  it  is  hopeless. 

The  results  of  simulations  for  the  designs  of  this  thesis  and  Capt 
Locken's  thesis  are  very  similar,  although  command  inputs  were  not 
identical.  Both  design  methods  provided  the  means  for  successfully 


compensating  the  aircraft. 

The  final  comment  in  the  comparison  of  the  two  methods  concerns 
compensation  bandwidth.  The  Porter  method  relies  on  straight  gains  in 
the  feedback  loops  to  achieve  desired  compensation.  Compensation  band¬ 
width  is,  therefore,  unknown.  In  the  Horowitz  approach,  the  bandwidth 
of  the  compensator  elements  is  a  prime  consideration  during  the  design 
process. 

Conclusions 

This  thesis  concludes  that  Professor  Horowitz'  quantitative  feed¬ 
back  technique  is  an  effective  tool  for  compensating  large  aircraft. 
The  aircraft  can  be  controlled  to  specifications  over  the  entire  range 
of  uncertainty  considered. 

The  sideslip  response  model  used  overestimates  the  speed  with  which 
the  aircraft  can  respond  to  a  bank  command.  Therefore,  the  bandwidth  of 
the  sideslip  loop  compensator  is  higher  than  necessary.  This  results  in 
sideslip  responses  which  far  exceed  the  specifications. 

Finally,  this  thesis  concludes  that  Professor  Horowitz'  method  does 
provide  inslgnt  into  the  overall  compensation  problem  throughout  the 
design  process.  An  example  of  this  is  the  ease  * n  analysing  the  reasons 
for  the  overdesign  in  the  sideslip  loop. 

Recommendations 

Quantitative  Feedback  and  its  application  to  flight  control  prob¬ 
lems  should  be  further  investigated.  A  more  complete  aircraft  model 
should  be  developed  to  include  bending  modes  and  coupling  between  the 
lateral  and  longitudinal  modes. 


Presently,  most  of  the  analysis  must  be  done  graphically,  by  hand. 
A  computer-aided-design  package  could  be  developed  to  greatly  enhance 
the  speed  and  accuracy  of  controller  design.  Ideally,  this  package 
would  incorporate  a  split-screen  representation  of  Bode,  Nichols,  and 
time  response  plots  which  change  as  the  user  changes  parameters  in  the 
transfer  function  of  the  model  or  compensator.  Computer-aided 
development  of  frequency  domain  bounds  based  on  response  models  and 
defined  range  of  uncertainty  would  also  be  useful. 
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APPENDIX  A 


Examples  of  Loop  Shaping 

Introduction 

The  fo lowing  loop  shaping  examples  are  provided,  courtesy  of  Dr. 
Horowitz.  All  Information  is  given  and  Nichols  plots  of  the  shaped  loop 
transmissions  are  provided. 


Example :  Shaping  of  a  nominal  loop  transmission  Lg(jw)  to  satisfy 
boundaries  B(w)  on  Nichols  Chart. 


Previous  notes  have  described  how  tolerances  on  the  closed-loop  system 
frequency  response  are  readily  translated  into  bounds  on  a  nominal  loop 
transmission  function  Lq(Ju).  In  Fig.  1,  for  example,  Lg(j2)  must  be  on 
or  above  the  curve  labelled  B(2),  etc.  is  the  "universal  high-frequency 
boundary"  applicable,  in  this  example,  to  i.e.  L^(ju)  (for  ui.40) 

must  be  contained  in  the  closed  curve  in  Fig.  1.  Additional  specifica¬ 
tion  is  eTa^>  where  •  is  excess  of  poles  over  zeros  of  L.(s).  Also,  L  is 
to  be  Type  1  (one  pole  at  the  origin).  We  proceed  to  describe  a  reasonable 
procedure  for  choosing  a  rational  function  Lq(s)  which  satisfies  the  above 
specifications. 

In  our  first  step,  we  try  to  find  the  B(u)  which  "dominates"  LQ(ju). 

E.g.  suppose  LQ1 (  j4)=0db/-135"  (point  A  in  Fig.  1).  But  at  w-l,  | LQ(j ) j 

needed  is  =27db.  In  order  to  decrease  |Lq|  from  27db  to  about  Odb  in  2 
2 

octaves  (4/1*2  ),  the  slope  of  |Lq(Juj)|  would  have  to  be,  on  the  average, 
about  -14db/octave,  Involving  yLg<-180a.  We  assume  "absolute"  stability 
is  required  here  for  Lg(ju)  with  a  margin  of  40°,  not  just  at  crossover 
("crossover"  is  defined  as  the  frequency  at  which  |Lg(joj) |*1).  Hence  B(l) 
dominates  Lg(ju),  at  least  more  than  B(4).  In  the  same  way  we  see  that 
B(l)  dominates  over  all  other  B(w)  in  Fig.  l.> 

The  B(u)  for  w<l  are  not  shorn  in  Fig.  1.  We  shall  assume  that  for 
uKl,  a  slope  of  -6  db/ octave  (with  27db  at  w“l,  i.e.  33  db  at  u*.5,  39db 
at  up. 25  etc.),  suffices.  We  can  tolerate  a-140°  for  u^l,  so  we  choose  a 
lag  corner  frequency  (symbol  lacf)  at  upl  (i.e.  pole  at  -1),  and  set 
|Lg(ju)|  (asymptotic)  at  30  db  (to  allow  for  the  -3db  correction).  Thus 
our  L.  is  so  far:  L_a31.6/s(s+l),  whose  phase  /Ln, (jw)  is  sketched  In 


'A'/ '.  .V.'.- 


*  • 

i 


>• 


S. 


yLg^(jq))  violates  the  -140®  bound  at  u>1.2,  so  a  lead  corner  frequency 

(symbol  lecf )  is  needed.  Where  should  it  be  located?  At  ui*5,  yLg^(j5)‘*-169® 

(see  Fig.  2),  so  29*  lead  is  needed;  but  we  know  that  later  there  will  be 

a  second  lacf .  so  allow  say  additional  15®  for  it  giving  15+29=45®  lead 

required  at  u>5;  which  is  achieved  by  a  lecf  at  w*5  i.e.  a  zero  at  -5. 

The  resulting  Lg2(s)*31.6(l+  j),  whose  phase  O‘0  is  sketched  in  Fig.  2. 

s(l+s)  i"~ 

In  the  Nichols  Chart  (N.C.)  we  are  (aw *10  or  so)  in  the  region  where  the 

maximum  phase  lag  allowed  is  135®  (i.e.  yLg(j(iQ  must  be  ^-135®).  Consider 

urlO,  with  present  /Lg^Om)  ■-112®,  so  135®-112®*23°  more  lag  is  allowed. 

But  this  lacf  will  be  followed  by  a  lecf,  so  allow  say  10°  for  it,  giving 

23°+10*»33o  more  lag  allowable.  This  locates  the  lacf  at  15.4  (tan33°*.65, 

and  10/. 65*15. 4) ,  so  we  set  the  lacf  at  i,i*15  (i.e.  pole  at  -15),  giving 

L03(s)-31.6(l+.2s)  /Lgi^)  is  sketched  in  Fig.  2. 
s(l+s)(l+|j) 

Looking  ahead  at  ui-40,  |Lg^(j40)j  =  -20db,  so  soon  Lg(ji  )  can  make  its 

asymptotic  left  turn  under  the  B.  boundary.  Our  plan  is  to  add  two  more 

"  "  n 

lecfs,  and  finally  two  complex  pole  pairs,  in  order  to  have  an  excess  e^ 

of  poles  over  zeros  of  4.  We  try  one  lecf  at  <i-.»40,  giving 

§ 

l04(«)-2imt±!  (1‘>‘  40)  .  /L04(ju)  is  sketched  in  Fig.  2. 

«(l+s)  (l+f^) 

We're  ready  now  for  the  last  lecf,  in  order  to  achieve  (an  asymptotic) 
horizontal  segment  for  lLg(Ju)|  ,  before  the  final  -24db/ectave  slope. 

(We  follow  Bode  in  this  respect,  a  good  master  to  follow.)  Where  should 
this  horizontal  segment  be  located?  The  bottom  of  (see  Fig.  1)  is  at 
-22.5db.  Allow  2db  margin,  3db  correction  due  to  the  last  lecf,  1.5db  for 
the  effect  of  the  lecf  at  '*'*40,  giving  a  total  of  -(22 . 5+2+3+1 . 5)*-29db. 
We'll  use  a  damping  factor  of  £*.6  for  the  2  complex  pole  pairs,  so  no 

A-4 


v 


correction  need  be  allowed  for  them.  Thus  the  final  break  for  |Lq(J(u)| 
asymptotic  is  to  be  at  -29db,  which  [L^(ju)|  achieves  at  w*60.  Hence, 
the  last  lecf  is  at  oi>60.  The  resulting  phase  due  to  and  the  lecf  at 
w*60,  is  -66°.  We  could  have  -180°  at  this  point,  but  we'll  allow  an  additional 
15°  margin  (a  matter  of  taste;  it  depends  on  the  problem  —  presence  of 
higher  order  modes,  etc).  This  means  100°  phase  lag  is  permitted,  50°  due 


to  each  complex  pole  pair  (180-66-15=100).  For  C-.6,  this  locates  them  at  LOO. 


Thus  Lq(s>* 


31.6(1+. 2s)  (1+  IqHI+Iq) 

*<i+s><i+fs>[i+Iuoi’+f5*>J 


Discussion 

l*Q(jw)  is  sketched  in  Fig.  1.  A  well-designed,  i.e.  "economical"  Lq(.iiu) 
is  close  to  its  boundary  B(«)  at  each  w.  The  vertical  line  -140°  is  the 


dominating  B(u>)  for  w<5  and  the  right  side  of  (line  -135°)  is  the  boundary 
effectively  for  5<u$30;  so  our  Lg(jui)  is  pretty  good  in  this  respect  since 
it  is  pretty  close  to  these  boundaries.  There  is  tradeoff  between  complexity 
of  Lq(s)  (number  of  its  poles  and  zeros)  and  its  final  cut-off  frequency,  now 
at  u*100.  There  is  some  phase  to  spare  between  ^(ju)  and  the  boundaries,  so 
use  of  more  poles  and  zeros  in  Lq(s)  would  permit  this  cut-off  frequency  to 
be  reduced  a  bit  below  100,  but  not  by  much.  On  the  ocher  hand,  if  we  want  to 


reduce  the  number  of  poles  and  zeros  of  Lq(s),  we  must  pay  the  price  in  a 
larger  cut-off  frequency.  We  could  economize,  significantly,  of  course, 
by  allowing  more  phase  lag  in  the  low  frequency  range.  If  -180°  was  permitted 
at  oj*l,  we  could  decrease  jLg(ju>)|  at  a  rate  of  12db/octave;  so  with  |Lq|* 

25db  at  w*l,  it  would  be  13db  at  01*2  (instead  of  the  present  18db).  Even 
with  no  more  saving,  this  5db  difference,  would  allow  a  cut-off  frequency 
at  about  70  Instead  of  100. 


Also,  Fig.  1  reveals  (immediately,  without  any  shaping  of  required) 
that  reduction  (i.e.  easing)  of  the  specifications  at  u*l  to  about  21db 
(instead  of  =26db),  would  have  the  same  effect  as  the  above.  One  can  check 
how  badly  the  specifications  are  compromised  by  such  easing.  The  design 
technique  is  thus  highly  "transparent"  in  revealing  the  trade-offs  between 
performance  tolerances,  complexity  of  the  compensation,  stability  margins, 
and  the  "cost  of  feedback"  in  bandwidth. 


DESIGN  EXAMPLE;  Unc«rtain  Plant  P-  k 

S(S+p) 

k:  [10*80] ,  p:  [-2.2]. 

QO,  ,  s} 

Note:  For  part  of  range  of  p,  plant  Is  unstable. 

W  -  t 


.pUrd-  te^rlc<T  -6 
O  p  l  d  U6 


Specifications  on 


[see  Fig 


k  a  I 


U&P. 


N.B.  At  high  a),  (A|t| )MAX  MUST  BE  >  |AP|MAX.  OK. 


Calculation  and  Construction  of  Plant  Templates  {P( 


The  templates  are  shovn  in  Fig.  B^.  It  suffices  to  calculate  (at  a 
fixed  ui)  several  values  of  l/J<s(jw+p) ,  i.e.  at  different  p  values.  This 


gives  the  bottoa  curve  of  the  template.  Then  extend  vertically  by 
18db(80/10-8;  201og8*18) . 


Step  2  Use  the  procedure  described  in  the  notes  (p.17  etc)  to  find  bounds 
on  Lg(jb))  in  order  to  satisfy  the  specifications  on  |T(ju)|.  A  nominal  plant 
must  be  chosen;  p»2,  k«10was  chosen  as  nominal  and  marked  heavily  in  Fig.  B0. 

a  — 

(It  helps  to  have  the  templates  on  transparent  paper  or  plastic.)  Already 
at  w»l,  the  bound  on  LqO)*  denoted  by  B(l)  is  determined  by  |L/(1+L)  (-63DB, 
rather  than  by  the  constraints  on  |T(Ju>) | .  And  this  is  so  for  w>l  also, 
which  is  not  typical  for  stable  plants  but  more  likely  for  plants  which  can 
be  open-loop  unstable  for  part  of  the  parameter  range. 


\ 


Shaping  of  Lg(jw) 

The  boundary  at  w*2  dominates,  i.c.  determines  the  level  of  |lq|  at 
ta)*2  to  be  slOdb.  The  boundary  at  ii>*5  determines  the  phase  there  to  be 
e-85°or  so.  This  necessitates  a  lead  corner  frequency  at  some  ^<5.  This 
was  chosen  to  be  at  1  and  lags  introduced  thereafter  such  that  at  u>«5  the 
phase  >  -85°;  they  were  chosen  at  5  and  8  (see  Figs.  B2>C).  A  lecf  is  then 
needed  to  be  followed  by  a  complex  pole  pair  with  £-0.6.  The  corner  below 
the  high-frequency  bound  can  be  turned  when  j L j  is  a  -25db.  If  we  try  a  lecf 
at  40,  then  phase  requirements  at  tu“20  force  the  final  cut-off  frequency  to  be 
well  beyond  the « value  at  which  |Lq|  *-25db.  A  lecf  at  w*25  gives  compat¬ 
ibility  of  phase  needs  at  m*20  and  turning  the  corner  when  jLgj  permits  it. 
This  gives  100  as  the  frequency  at  which  the  final  complex  pole  pair  can  be 
inserted.  This  gives  Lg(s)«  4(1+s)(1+|j) 

»<l+.2s)a+  |)(1+  +  (foo)2’ 

*GP.  ■  »  giving  G(s),  —  see  Figs.  B,,,  C. 

8 (s+2) 

Finding  F(s) 

The  proper  G  guarantees  thati|T|  does  not  exceed  those  allowed.  The 

next  step  is  to  find  the  range  of  jL(ji4i)/(l+L(j(ii)  |.  Place  the  template  of 

the  plant  at  w* 5  (for  example)  on  the  point  Lg(j5)  in  the  Nichols  Chart, 

i.e.  on  5db  / —82?  and  it  is  seen  that  |L/(1+L)[  max»  2.4db,  min  -1.7db. 

But  the  specifications  require  |T  j»  |FI./ (1+L)  |e  [-35,  -17]  db.  Therefore,  it 

is  required  that:  -33. 3db< | F|<  — 19 . 4db.  In  this  wav  we  obtain  the  bounds 

on|F(ju)|  shown  in  Fig.  D.  It  is  easy  to  find  an  F(s)  which  satisfies 

these  bounds,  e.g.  F(s)  ■  _ 1 _ 

(s+nu+^m+fg)  . 
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APPENDIX  B 


Introduction 

This  appendix  contains  the  output  from  the  computer-aided-design 
program,  TOTAL,  of  the  response  modeling  process.  Included  are  the 
transfer  functions,  the  time  responses,  and  the  frequency  responses  for 
the  bank  angle  and  the  sideslip  angle  models. 


CLOSED-LOOP  transfer  function 
BANK  ANGLE  UPPER  BOUND  (Bll) 


CLK*  <  CLNK/CLDK  )*  5.800 

CLTF<S>  NUMERATOR 

I  CLNPOLY< I )  CLZERO(I) 

1  <  5.800  >  CLNK*  5.888 


CLTF<S>  DENOMINATOR 


I 

CLDPOLY ( I > 

CLPOLE< I > 

1 

<  1.008  )S»*  2 

•:  -.5000  >  +  J<  0. 

> 

2 

<  10.50  >S**  1 

<  -10.00  >  +  J-:  0. 

> 

3 

■:  5.000  > 

CLDK*  1.000 

OPTION  >37 

CONTINUOUS  TIME  RESPONSE  FOR  CLTF'S) 
WITH  STEP  INPUT  OF  STRENGTH  =  30. 


RISE  TIME: 
DUPLICATION  TIME: 
PEAK  TIME: 
SETTLING  TIME: 
PEAK  VALUE: 

FINAL  VALUE: 


TR*  4.39982 
TD*  VERY  LARGE 
TP*  VERY  LARGE 
TS*  7.92643 
MP*  38.0008 
FV~  38.0000 


OPTION  >31 

CONTINUOUS  TIME  RESPONSE  FOR  CLTF<S> 

WITH  STEP  INPUT  OF  STRENGTH  *  38 . 

ENTER  INITIAL  TIME,  FINAL  TIME,  DELTA  TIME  >0,15,1 


0.  0. 


1  .0000000 

10.844472 

2.0000008 

18.382754 

3.0000000 

22.953784 

4.0000000 

25.726254 

5.0000000 

27.407842 

6.0080000 

28.427777 

7.0000000 

29.046398 

8.0000000 

29.421411 

9.0000000 

29.649190 

10.000000 

29.787223 

11.000000 

29.870944 

12.000000 

29.921724 

13.000000 

29.952523 

14.000800 

29.971294 

15.000000 

29.982534 

B-3 


OPTION  >32 

CONTINUOUS  TIME  RESPONSE  FOR  CLTF<S)  WITH  STEP  INPUT  OF  STRENGTH  »  30 
REGION  OF  CALCULATION:  T-  0.  TO  T«  10.00 

F<T>*  0.  TO  F(T)«  30.00 

0 

#44444444444444444444444444444444444444444444444444444444F(T) 

4  * 

4  »  - 

4  # 

4 


*  - 


CLOSED-LOOP  FREQUENCY  .RESPONSE  USING  OPTION  54 
MAGNITUDE  <**#>  AXIS  RANGE:  -70.0  TO 

FREQUENCY  < RAD/ SEC) 

.112  ♦  •  -  -  * 

.126  +  -  -  * 

.141  ♦  -  -  -  * 

.138  +  -  -  -  * 

.173  ♦  -  -  -  * 

.200  ♦ - *— 

.224  +  -  -  # 

.231  +  -  -  -  * 

.282  ♦  -  -  -  * 

.316  ♦  -  -  -  * 

.353  * - * — 

. 398  +  -  -  -  * 

.447  +  -  -  -  * 

.301  ♦  -  -  -  * 

.362  + - - - - - - — * 

.631  +  -  -  -  * 

.709  +  -  -  -  * 

.794  + - * - 

.891  +  -  -  -  * 

1.00  - - * - 

1.12  ♦  -  -  -* 

1.26  ♦  -  -  -* 

1.41  ♦  -  -  * 

1.38  *  -  -  * 

I. 78  +  -  -  *- 

2.08  ♦ - * - 

2.24  ♦  -  *  - 

2.31  +  -  *  - 

2.82  ♦  -  *  - 

3.16  +  -  *  - 

3.33  + - * - 

3.98  ♦  -  * 

4.47  +  -  -  *  - 

3.01  +  -  -  * 

3.62  + - * - 

6.31  +  -  -  # 

7.08  ♦  -  -  * 

7.94  + - * - 

8.91  +  -  -* 

10.0  ♦ - * - 

II. 2  *  -  * 

12.6  +-•»-- 

14.1  +  *  - 

13.8  ♦  *  - 

17.8  +  # 

20.0  + - * - 

22.4  ♦  * 

23.1  +  * 


DECIBELS 


,•  ,•  *. *  *,*  ■  -  *  -  ■„«  •  ‘.•"j  '  -  » •  .  • .  ,’•*  »*■ 


CLOSED-LOOP  TRANSFER  FUNCTION 
BANK  ANGLE  LOWER  BOUND  <A11> 


<  CLNK/CLDK  >= 


CLNPOLY< I ) 
5.068  ) 


CLTF<S)  NUMERATOR 


CLZERO< I ) 
CLMK*  5.00 


CLTF<S)  DENOMINATOR 


CLDPOLY<I) 

I . 080  >S**  3 

II . 50  >S»#  2 

15.50  >S«*  1 

5.000  ) 


CLPOLE<I> 

-.5808  >  +  J<  0 

-1  .008  )  +  J(  0 

-10.00  >  +  J<  0 

CLDK*  1.00 


OPTION  >37 

CONTINUOUS  TIME  RESPONSE  FOR  CLTF<S) 
WITH  STEP  INPUT  OF  STRENGTH  *  38. 

RISE  TIME:  TR»  5.18511 

DUPLICATION  TIME:  TD*  VERY  LARGE 

PEAK  TIME:  TP»  VERY  LARGE 

SETTLING  TIME:  TS»  9.30233 

PEAK  VALUE:  M P*  30.0000 

FINAL  VALUE:  FV»  30.8000 


OPTION  >31 

CONTINUOUS  TIME  RESPONSE  FOR  CLTF(S) 

WITH  STEP  INPUT  OF  STRENGTH  *  38. 

ENTER  INITIAL  TIME,  FINAL  TIME,  DELTA  TIME  >8,15,1 


0. 

I . 3080000 
2.0000000 
3.0000000 
4.0000000 
5.0000000 
6.0000000 
7.0000000 
8.0000000 
9.0000000 
10.000000 

II. 000000 
12.000000 
13.000000 
14.000000 
15.000000 


0. 

3.9554485 
11 .276685 
17.567138 
22.063038 
25.048233 
26.938179 
28.123193 
28.854485 
29.382493 
29.575959 
29.742445 
29.843652 
29.905121 
29.942435 
29.965079 


CLOSED-LOOP  TRANSFER  FUNCTION 
SIDESLIP  ANGLE  UPPER  BOUND  <821) 


CLK*  <  CLNK/CLDK  >•  2.000 


1 
1 

2 
3 


CLTF<  S)  NUMERATOR 


CLNPOLY< I ) 

<  2.090  >S**  2 

<  1.400  >S»*  1 

<  0.  ) 


CL2ERO< I ) 

0.  )  +  J<  0. 

-.7000  )  +  J<  0. 

CLNK*  2.000 


CLTF<  S)  DENOMINATOR 


I  CLDPOLY<I> 

1  <  1.000  >S**  3 

2  <  24.00  )S»*  2 

3  <  83.00  >8**  1 

4  <  60.80  ) 


CLPOLE<I) 

<  -1  .000  )  +  J<  0.- 

<  -3.800  )  +  J<  0. 

<  -28.00  )  +  J<  0. 

CLDKa  1.000 


OPTION  >37 

CONTINUOUS  TIME 
WITH  STEP  INPUT 

RISE  TIME: 
DUPLICATION  TIME: 
PEAK  TIME: 
SETTLING  TIME: 
PEAK  VALUE: 

FINAL  VALUE: 


RESPONSE  FOR  CLTF<S> 
OF  STRENGTH  *  38. 

TR-  UNDEFINED 
TD*  UNDEFINED 
TP*  .107202 

TS*  VERY  LARGE 
MP*  2.09694 

FV*  8. 


OPTION  >31 

CONTINUOUS  TIME  RESPONSE  FOR  CLTF<S> 

WITH  STEP  INPUT  OF  STRENGTH  »  30. 

ENTER  INITIAL  TIME,  FINAL  TIME,  DELTA  TIME  >0,15,1 
0.  0. 


1  .0008000 

.27818235E-01 

2.0000000 

-.54045369E-81 

3.0000000 

-.23882450E-01 

4.0000000 

-.36508907E-02 

5.0000000 

-.3190417SE-02 

6.0000000 

•  -. 1 1748840E-02 

7.0000000 

-.43194101E-03 

8.0000000 

-. 15890320E-03 

9.0000000 

-.58457268E-04 

10.000000 

-.21505230E-04 

11.000000 

-.791 13319E-05 

12.000000 

-.29104164E-05 

13.000000 

-.  10706324E-05 

14.000000 

-.39388282E-06 

15.000000 

-.14498] 10E-96 

8-9 


/  m'  •  -  « 


«P 


CONTINUOUS  TIME  RESPONSE 
REGION  OF  CALCULATION:  T= 

FCT)= 

8 


FOR  CLTF<S)  WITH  STEP  INPUT  OF  STRENGTH  =  30. 
0.  TO  T=  10.00 

-.4000  TO  F(T)=  2.000 
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APPENDIX  C 


Trans-fer  Functions 


Introduct i on 

This  appendix  contains  the  trans-fer  -functions  o-f  the  e-f- fectiue 
plants  used  in  the  design  o-f  the  bank  angle  loop  and  the  trans-fer 
■functions  o-f  the  compensated  system  used  during  the  simulations  o-f  the 
design . 


E-f-fective  a's 

Cruise  1 


<111 


* 


_ 9.7278  <t  *  22.35)  (s  ♦  15 .66) _ 

<s  +  9.5295)  <s  +  23.27)  <s  +  16.52)  (s  -  9.91572)  <C-1 > 


Cruise  2 


<111 


* 


_ 9.7983  <s  +  8.3731)  (s  +  39.3) _ 

<s  +  9.785  +  j  9.941)  <s  +  28.71)  <s  -  9.932) 


Power  Approach 


<»11 


# 


1.434  (s  +  8. 996)  <  s  +  1.37)  (s  *  4.12) 
s  <s  +  1.28  +  j  9.23)  <s  +  3.42)  <s  -  9.94) 


■:C-3) 


Simulation  Trans-fer  Functions 
Cruise  1 


Denominator 


< s+9 . 8+ j 9 . 3)<s+ 1  . 2)  < s+8 . 35)  < s+9 .7)  <s+3. 7+j 4.3)-: 
x  <  s+8 . 4+ j  5)  <  s+23 . 3) ( s+ 16 .5) <  3  +  96 . 6+j  1 1  9 .5)  s+ 


s+ 5 . 3+..:  7.7) 
125+jl“6> 


Numerators 


Bank  Angle  Response 

26??  <5+1  .2Xs+l  Xs+8 .86X s+S. 7+ j5)  <s+22.4Xs+15.7) 
x  <s+96.7+jll9.6Xs+125+jl76> 

Control  Wheel  De-flection 

3788  <  s+ 1  X s+ 1 .2)<5+8.87)<s+8 .52)<s-0.016Xs+9.3+j5.2) 
x  < 5+22.3)  <  s+14 .5)  <  5+96 . 8+,  j  1 19.5) < S+125+.  j  176) 

Sideslip  Response 

85.9  s< s+l X 5+8. 86) < 5+2) < s+2. 4) < s+20) <5+11 4+ j 152) 2 
x  <s+23.3) 

Rudder  De-flection 

-5.9x1013  <s+l )3<s+0 .06) <5+2.4) < s+5) <5+22.4) 

Cruise  2 

Denominator 

<5+0. 9)  <5+0  .7X5+1  .lXs+0.96+j0.65Xs+3.6+j4.3) 
x  <  s+5. 4+ j7. 7Xs+6. 5X  s+14. T+j  14. 7X5+30. 3) 
x  < s+95 . 5+j 1 1 5) < s+ 1 25 . 3+ j 177.7) 

Numerators  * 

Bank  Angle  Response 

2630  <5+0 .9X s+l  X 5+1 .15X5+6.4)  < s+l 4+ j  15.2)  <5+38 .3) 
x  <  s+95 . 7+j 1 1 5 . 2)<  5+ 1 25 . 3+ j 1 77 . 7) 


Control  Wheel  Deflection 


3708  <  s+1  Xs+  .89Xs+0  .?)< s-0 .83) < s+1 .12X3+6.5) 
x  <s+i4+jl4.7Xs+30.3Xs+95.5+jll5Xs+125.3+jl77.7)  <C-1 1  > 


Si desl i p  Response 


-131.6  s<  3+0.095X5-1 .  6X  s+1  X  s+2X  s+30 . 3X  s+20 ) 
x  <s+ll  4+ j  132)  <  5+114+^152)  <C-!2) 


Rudder  Deflection 


1.04X10*4  <s+0.095Xs+l)3<s-l  .6Xs+5Xs+30 . 3)  <C-13> 


Power  Approach 


Denominator 


<s+0 .8+,j0 .03) <s+0 ,7)< s+0  .*006) <s+l  ,3Xs+2.3+j2.5) 
x  <s+T.5+j3.7Xs+4.3+ j8.82Xs+5.6+ j8.1Xs+49.1> 
x  <s+82.6+jl87.6Xs+T28.5+jl81 .6)  <C-14) 


Numerators 


Bank  Angle  Response 


5318  < s+1 ) < s+0 .8) < s+1 .4) <s+0 .006)<s+2.5+j2.2) 
x  <s+4.1Xs+49.5Xs+82.6+jl87.6Xs+128.3+jl81  .6)  <C-13) 


Control  Wheel  Deflection 


3788  s<s+l  .l)<s+l)<s+0.8Xs+l  .6X5-0.04X3+2.3+ J2.3) 
x  < s  +  4 . 1  X s+49 . 1 X s+82 . 6+,.j  107.6)  < s  +  1 28 .5+,j  131 .6)  <C-!6  > 


Sidesl ip  Response 

-16.9  s2'  s+8 . 09X  s+ 1  >  <  s+2)  '  s+4 . 1 )  ■'  s+2?  '  •'  s-63^ 
x  <  s+ 1 14+_j  152)  <  5*  1 1 4+^j  1 32) 
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